DUDLEY K!iOX LI s:\MH 
NAVAL PCSTC 'OUAin SCL'CL 
MONTEREY. C/'LIF0R.;ta rS''43 



NAVAL POSTGRADUATE SCHOOL 

Monterey, California 




THESIS 

ALGORITHMS AND HEURISTICS FOR 
TIME -WINDOW-CONSTRAINED 
TRAVELING SALESMAN PROBLEMS 

by 

Chun, Bock Jin 
and 

Lee, Sang Heon 
September 1985 

Thesis Advisor: Richard E. Rosenthal 

Approved for public release; distribution is unlimited. 

T226058 





^ »• 















security classification of this page rw»#fi Dmtm Bntmfd) 



REPORT DOCUMENTATION PAGE 


READ INSTRUCTIONS 
BEFORE COMPLETING FORM 


1. REPORT NUMBER 


2. GOVT ACCESSION NO. 


3. RECIPIENT'S CATALOG NUMBER 


4. TITLE 5u6f/f]©; 

Algorithms and Heuristics for Time-Window-Con- 
strained Traveling Salesman Problems 


S. TYPE OF REPORT 4 PERIOD COVERED 

Master's Thesis; 

September 1985 


6. PERFORMING ORG. REPORT NUMBER 


7. AUTHORfs; 

Chun, Bock Jin 
Lee, Sang Heon 


8. CONTRACT OR GRANT NUMBERr*; 


9. PERFORMING ORGANIZATION NAME AND ADDRESS 

Naval Postgraduate School 
Monterey, California 939^3-5100 


10. PROGRAM element. PROJECT, TASK 
AREA 4 WORK UNIT NUMBERS 


11. CONTROLLING OFFICE NAME AND ADDRESS 

Naval Postgraduate School 
Monterey, California 939^3-5100 


12. REPORT DATE 

September 1985 


13. NUMBER OF PAGES 

103 


14. MONITORING AGENCY NAME 4 ADDRESS(lf different from Controlling Office) 


IS. SECURITY CLASS, (of this report) 


1S«. DECLASSI FICATION/ DOWNGRADING I 

schedule 



16. distribution statement (of this Ropoet) 



Approved for public release; distribution is unlimited. 



17. DISTRIBUTION STATEMENT (of the abstrect antaeed In Block 20, If different from Report) 



18. SUPPLEMENTARY NOTES 



19. KEY WORDS (Continue on reverse side If necessary and Identify by block number) 

Heuristic, Algorithm, Time Window, Hard Time Window, Soft Time Window, 
Slack, Branch and Bound, Nearest Neighbor, Penalty Cost, Traveling 
Salesman Problem, State-Space Relaxation 



20, ABSTRACT (Continue on reverse elde If necessary and Identify by block number) This thesiS report S On m^^th — 

ods for solving traveling salesman problems with time-window constraints. Two 
types of time windows are considered: hard time windows, which are inviolable, 
and soft time windows, which are violable at a cost. For both cases, we devel- 
op several heuristic procedures, including some that are based on Stewart's 
[Ref. 6] effective heuristics for the traveling salesman problem without time- 
window constraints. In addition, we develop exact algorithms for each case, 
which are based on the state-space relaxation dynamic programming method of 
Christof ides , Mingozzi, and Toth iRef. 5] . Computational experience is re- 
ported for all the heuristics and algorithms we develop. 



DD 1 jan*^73 1473 edition OF 1 NOV 65 IS OBSOLETE ^ 
S N 0102- LF- 014- 6601 



SECURITY CLASSIFICATION OF THIS PAGE (When Data Entered) 



Approved for public release; distribution is unlimited. 



Algorithms and Heuristics for . 
Time -Window-Constrained 
Traveling Salesman Problems 

by 

.Chun, Bock Jin 

Major, Republic of Korea Air Force 
B.S., Korea Air Force Academy, 1976 

and 

Lee , Sang Heon 

Major, Republic of Korea Army 
B.S., Korea Military Academy, 1977 



Submitted in partial fulfillment of the 
requirements for the degree of 

MASTER OF SCIENCE IN OPERATIONS RESEARCH 

from the 

NAVAL POSTGRADUATE SCHOOL 
September 1985 



ABSIB&CT 



This thesis reports on methods for solving traveling 
salesman problems with time-window constraints. Two types of 
time windows are considered: hard time windows, which are 

inviolable, and soft time windows, which are violable at a 
cost. For both cases, we develop several heuristic proce- 
dures, including some that are based on Stewart’s [Ref-6] 
effectxve heuristics for the traveling salesman problem 
without time-window constraints. In addition, we develop 
exact algorithms for each case, which are based on the 
state-space relaxation dynamic programming metnod of 
Christo fides, Mingozzi, and Toth [Eef-5]. Computational 
experience is reported for all the heuristics and algorithms 
we develop. 
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I- INIflODOCTI ON 



A- OVEBVIEH 

Consider a traveling salesman having to visit n cities 
or customers. He starts from a depot and needs to visit each 
of the other n-1 cities only once and then return to the 
depot. The ccst of traveling between any pair of cities 
(expressed in terms of distance, time or cost, etc), say 
from city i tc j, is given as crj in a cost matrix C. The 
problem is to design a tour through the n cities that would 
minimize the total cost of the tour. This is known as the 
Traveling Salesman Problem which is a well-known classical 
operations research problem. 

The TSP is called Euclidean when the cities that must be 
visited all lie on the same plane and the cost of traveling 
between any pair of cities is the Euclidean distance between 
them. 

The TSP is an NP-complete problem [ Eef . 1, 2]. All known 
exact solution methods have a rate of growth of the computa- 
tion time which is exponential in c. On the other hand, 
heuristic solution methods have a rate of growth of the 
computation time which is a low order polynomial in n and 
have been experimentally observed to perform well. For this 
reason, there has teen an extensive amount of research 
directed at TSP heuristics. 

In this thesis we consider adding time window 
constraints to the TSP. That is, if t; is the time tna t the 
salesman visits city i, then t;,mast satisfy 1^ ^ < u- , 

where li and Uj^ are the specified lower and upper bounds of 
a time window. This problem is not as well studied as the 
unconstrainted TSP, but there have teen a few approaches 
used on the problem. 
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Psaraftis [Ref. 3] has presented a dynamic programming 
model and solution procedure for two dial-a-ride problems, 
which are similar to time-window constrained TSPs- Baker 
[Ref. 4] has presented an exact algorithm using a branch and 
bound procedure which is effective for very small n. 
Christofides et al. [^Ref. 5] have presented a dynamic 

programming state space relaxation procedure to compute 
bounding information within a branch and bound algorithm. 

The objective of this study is to develop exact and 
heuristic algorithms which will provide an optimal or near 
optimal tour that visits each city in its given time 
window. He are given a depot location, a set of x,y 
co-ordinates for n cities and a set of time windows . 

a common application of the TSP is in vehicle routing 
problems. A set of customer orders must be partitioned 
among several vehicles. Given a partition, the problem then 
decomposes into one TSP for each vehicle. Because of this 
prospective application and in deference to the difficulty 
of large TSPs (with or without time constraints), we confine 
our research and computation to small-scaled problems (less 
than 30 nodes). 

We consider two different kinds of time window 
constraints: hard time windows and soft time windows. Hard 

windows cannot be violated. Soft windows can oe violated, 
but a penalty cost must be paid when they are. The penalties 
can be defined individually for each customer, and they can 
differ for early and late arrivals. Generally, the penalty 
for arriving before the lower time window bound is much less 
than the penalty for arriving after the upper bound. In 
Chapter III, we present the hard time window approach and in 
Chapter IV, we present the soft time window approach. 

He developed several Fortran programs for solving the 
TSP and time-constrained TSP. For the TSP, we use Stewart's 
[Ref. 6] recent heuristics, CCCO and CCAO. For the time 
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constrained ISP problems, we develop some new heuristics, 
seme of which are modification of Stewart’s heuristics for 
the unconstrained problem. We also developed exact algo- 
rithms for both hard and soft windows using Christofides et 
al. *s [Ref. 5] method of state-space-relaxat ion dynamic 
programming and branch and bound. This is described in 
Chapters III and IV. 

Finally, we describe a hybrid of the heuristic and exact 
programs. The hybrid uses the overall structure of the 
exact program, bat the upper bounds are obtained with the 
heuristic. Ihis is discussed in Chapters III and IV. 

B. THE TSAVEIIllG SALESHIB PBOBLBH 

A tour is a chain which passes through all the n nodes 

and in wnich the first and the last nodes coincide. A tour 

is also known as a Hamiltonian cycle. 

Let a tour be denoted by t = (i , i ,..., i , i ) and 

12 n 1 

the cost of this tour be 
n-1 

C (t) = c + c . 

i i i i 

j = l j j+1 n 1 

Here <i,i,...,i)is a permutation of the inteaers 
1 2 n 

from 1 to n, giving the order in whica the cities are 
visited . 

The Traveling Salesman Problem can be defined as 
follows. Given a graph G = { N',A } composed of a set of 
nodes N, a set of arcs A connecting these nodes, and a cost 
(distance ) c^j associated with each arc (i, j) in A. The 
ISP is the problem of finding the minimum cost tour of the 
nodes in N. The following mathematical formulation of the 
TSP is from Stewart [Ref. 6]. 
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ij 

y . > 0 for all (i, j) (1,7) 
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where 

X = / 1 arc (i, j) is on the tour 

i j ( 

\ 0 otherwise 

y are continuous variables that force the final 
ij 

soluticn to be on the tour 

( i.e. include every node in the same route ) 
u is a number > n- 1 , and 

n is the number of nodes in the set N. 



The constraints (1.2) and (1.3) ensure that each node 
will be visited exactly once, while constraints (1.4) and 
(1.5) force the final solution to be a single tour that 
starts and ends at node 1 (depot) . This formulation is not 
directly used in our TSP programs, but is of interest in a 
general discussion of the problem- 
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C. TSP WITH TIHE WIHDCW COHSTEAINTS 



The time-ccnstrained Traveling Salesman Problem is a 
variation of the TSP that includes time window constraints 
on the time to visit some of the cities. The hard time- 

constrained TSP is to find the minimum cost tour subject to 

visiting each city within its time window. 

‘For the time-constrained TSP model, we define a contin- 
uous nonnegative variable, t^ , to be the time that the 
salesman visits city i. Since the salesman must return to 
city 1 (depot) at the end of the tour, the formulation 

includes an additional variable, t„+| , the total time 

required to complete the tour. 

Me assume that a complete, symmetric, nonnegative 
distance matrix, |Cij |, is known and that time is a scalar 
transformation of distance. Thus time and distance may be 
used interchangeably. 

The following mathematical formulation of the TSP with 
time constraints is from Baker [2ef. 4]. 



MIN t 

n+1 


t 

1 






(1.8) 


S.I t 

i 


t > 
1 


c 

1i 




(1.9) 


1 t, - 

X 


t j > 
j 


c 

ij 


j = 3^4, 


(1.10) 








2 < i < j 




t - 
n + 1 


t > 

i 


c 

1i 


i — 2/3, 


(1-11) 




t > 

i 


0 


i = 1/2,... ,n+1 


(1.12) 


VI 


VI 




X ^ 2,3,. ../O 


(1.13) 



i i i 

where t = the time that the salesman visits city i 
i 

|x| = the absolute value of x 

c = the shortest time required to travel from 

ij 

city i to city j 
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1 = the lower bound on the time window for the 

i 

salesman to visit city i 

by assumption all 1 > 0 

i 

u = the upper bound on the time window for the 
i 

salesman to visit city i 

u > 1 , for all i 

i i 

The constraints (1-9) through (1.12) ensure a nonnega- 

tive arrival time at city i, t- , be obtained for each city 

(node 2 through node n) on the tour. The constraint (1.9) 

guarantees that t; , the time that the salesman leaves the 

node 1 (depot) will be the smallest t value. The absolute 

value constraint (1.10) ensures that the arrival times of 

any two city differ by amount of time sufficient to allow 

the salesman to travel between the two cities. The 

constraint (1.11) guarantees that t^+i , the time the salesman 

returns to the depot, will be the largest t value- The 

inequalities (1.12) and (1.13) are nonnegativity and the 

time window constraints respectively- 

Unf ortunately. Baker's proposed model for the time- 

constrained TSP is very difficult to solve, because 

constraint (1.10) is nonconvex. Therefore, we will not use 

this formulation directly in our program. 

Figure 1.1 illustrates the nonconvexity of constraints 

(1-10) for one i,j pair, the example | t^ - t^ | > 5. The 

feasible region for this constraint is the union of two 

disjoint sets. Taken all together, constraints (1.10) 
m 

define 2 disjoint sets sets where m = (n-1) (n-2)/2, 

which are very difficult to work with. 

We can see that the time-constrained TSP is very 
different from ISP, even in formulation. 
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Figure 1.1 Exaaple of Nonconvexity of (1.10) 
in Two Dimensions. 
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II. HEO B IST IC TSP SOLOTION 



A- OVEB?IEW 

Many heuristic procedures have been developed for 
solving TSP. Our purposes in this Chapter are to examine 
some of the well-known heuristics, to review Stewart's [Ref. 
6] recent heuristic, and to compare these approximate tech- 
niques on the basis of efficiency and accuracy on a small 
number of exaiples. 

In general, heuristic procedures are categorized by 
three broad classes: tour construction procedures, tour 

improvement procedures, and composite procedures [Ref. 7]. 
Tour construction procedures start with a single node and 
successively add nodes till a tour is built. Tour improve- 

ment procedures attempt to find a better tour given an 
initial tour. Composite procedures construct a starting tour 
from one of the tour construction procedures and then try to 
find a tetter tour using one cr more of the tour improvement 
proced ures. 

1- Tour Con struction Procedures 

There are many methods available for constructing an 
initial tour. Procedures which have been generally used are 
given below. 



a. Nearest Neighbor { Rosenkrantz et al. [Ref. 3] ) 

Step 1. Start with any node as the beginning of 
a subtour. 

Step 2. Find the node closest to the last node 
added to the subtour. Add this node to 
the current subtour. 
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step 3. Repeat step 2 until all nodes are contained 
in the tour. Then, join the first and last 
node. 



b. Clarke and Wright Savings ( Clarke and Wright 

[fief- 9] ) 

Step 1. Select any node as the central depot 
which we denote as node 1. 

Step 2. Compute savings s = c + c - c 

ij 1i 1j ij 

for i,j = 2,3,...,n. i?^j 

Step 3. Order the savings from largest to smallest. 

Step 4. Starting with the largest savings on the 
list, subtours are assembled such that the 
next node added has the largest remaining 
savings - provided that a constraint is not 
violated. Once a pair of nodes i and j have 
been linked, they remained linked . 

Repeat until all nodes have been assigned. 

Here, tne quantity s, _ is the amount of travel 
saved if node j is visited direc?.ly after i, as opposed to 
having separate trips from the depot to nodes i and j. 
Figure 2.1 demonstrates the procedure for two nodes i and j. 

c. Insertion procedures ( Rosenkrantz et al. [Ref. 

8 ] ) 

An insertion algorithm constructs a feasible 
tour by successively adding one node to an existing subtour. 
This procedure takes a subtour of k nodes at iteration k ana 
attempts to determine which node not in the subtour should 
join the subtour next {the Selection step) . Then it deter- 
mines where in the subtour it should be inserted (Insertion 
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Figure 2.1 Concept of the Clarke - Bright 
Savings Beuristic. 



step). Stewart [ Eef . 6] presented the following general 

algorithmic structure. 



Step 1. (Initial Suttour) 

Octain a TSP tour for a subset of the nodes 
N'c N in G. 

Step 2. (Selection Step) 

Find a node k € N-N * to be added to the 
existing subtour. 

Step 3. (Insertion Step) 

Choose an arc(i,j) in the subtcur on N'. 
Insert node k between i and j and add 



k to N' . 

Step 4 . If N = N', then stop- 

(We have a Hamiltonian cycle). 
Otherwise, return to step 2. 
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There are many variations on this algorithmic 
structure depending on the procedures chosen for executing 
steps 1 ,2 and 3 . 

Wiorkowski and McElvain [Ref. 10], Or [fief. 11], 
Stewart [Ref. 12] and Norhack and Love [Ref- 13] all present 
insertion algorithms that use the convex hull of the set of 
nodes N for the initial subset N*. Nemhauser and Hardgrave 
[Ref- 14] have shown that there exists an optimal tour for 
every Euclidean TSP in which the relative order of the nodes 
on the boundary of the convex hull is preserved. This means 
that the optimal tour visits nodes on the boundary of the 
convex bull in the same order as if the boundary itself were 
followed. 

Further justif ication for the use of the convex 
hull for the initial sub tour is shown empirically by 
Stewart's [Ref. 6] computational experiment. He compared 
several insertion heuristics both with and without the 
convex hull as the starting solution. The results show that 
all the insertion algorithms are improved by the use of the 
convex hull. Some are improved substantially, others only 
moderat ely. 



Many criteria have been suggested for the selec- 
tion of the node to be inserted in an insertion procedure. 

( 1) N earest Neighbo r ( Rosenkrantz et al. [ Ref. 
8] ). Choose the node k that is nearest a node in the 
current tour. I. e. , find k = argmin c s.t. j€K-N*, 

(2) Cheapest I n ser tion ( Ro s en krantz et al. 
[Ref. J] ) . Choose the node k that may be inserted at 
minimal increased cost. I.e., find 

k = arg min c +c -c s. t. m£N-N*, i,j£N'. 

m im mj ij 
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(3) Farthest Insert ion ( Bos en kr antz et al. 
f Bef . 8]). Choose the node k that is farthest from a node 
current subtour. I.e., find k = argmax c s.t. j 6 N-N', 



r Ref. 8] ) 



Arbitrary Insert ion { Rosen k rant z et al. 

Choose node k randomly from among M-N*. 

(5) Ratio I nser tion ( Stewart [Ref. J2] ) . 

Choose the node k such that the proportional increase in 

cost is minimal. I.e., find k = argmin (c + c ) / c 

m im mj ii 

s.t. m € N-N • , i, j € N* . 



(6) Pe rpendicul ar Dist^ce ( Wiork owsk i and 
McSlv ai n [Ref. JO] ) . - Choose the node k that is closest to 
an arc in the current subtour. 

(7) Rati o Ti mes Distance ( Or [Ref. JJ] ). 

Choose the node k such that the product of ratio and 

distance is minimized. I. e. , find 

k = argmin ((c +c )/c )x(c +c -c ) 
m im mj ij im mj ij 

s.t. m € M^N* , i,j e N*. 

(fi) Greatest An^le ( No rba ck and Love [Eef. 

13 1 ). Choose the node k and arc i, j such that the angle 

formed by the two arcs (i,k) and (k,j) is a maximum. I.e., 

find k = argmax angle{ arc(i,m), arc(m,j) } s.t. m£N-N*, 
m 

i/j 6 N* . 



The insertion criteria that have been used 
fall into two categories. [Ref. 6] 



1. Cheapest Insertion 

Insert the node k6 N-N' between those two 

connected nodes i,j € N' that minimize the guantity 

c + c - c . 
ik kj ij 
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2. Identical Insertion and Selection 



Do selection and insertion in the same 

step. 

2. Tou r Improvement Procedures 

The best known procedures of this type for the TSP 
are the branch exchange heuristics [Ref- 7]. These branch 
exchange heuristics work as follows. 

Step 1. Find an initial tour. This tour may be chosen 
randomly from the set of all possible tours, 
cr it may be generated by one of the tour 
building procedures above. 

Step 2. Improve the tour using one of the branch 
exchange heuristics. 

Step 3. Continue step 2, until no additional 
improvement can be made. 

For a given k, we define a k-change of a tour as 
consisting of the deletion of k branches in a tour and their 
replacement by k other branches to form a new tour. A tour 
is k-opt if it is not possible to improve the tour via a 
k-change. In general the larger the value of k, the more 
likely it is that a k-opt solution will be optimal. 
Unfortunately, the number of operations necessary to test 
all k exchange is proportional to n , where n is the number 
of nodes in the TSP. Due to this complexity, values of k = 2 
and k = 3 are most commonly used [ Ref. 7]. The 2-opt and 

3-opt heuristics were introduced by Lin [Ref. 15] and tna 
k-opt procedure, for k>3 was presented by Lin and Kernighan 
[ Ref. 16]. 

Or [Ref. 11 ] has designed a modified 3-opt that 

considers only a small percentage of 3-branch exchanges. 
This modified 3-opt called Oropt by Stewart [Ref. 6] 



21 



considers onlj those branch exchanges which are composed of 
a string of one, two, or three adjacent nodes being inserted 
between two ether nodes in the current tour. 3y limiting the 
number of exchanges that are considered in this way, Oropt 
requires many fewer calculations than a full 3-opt. 

Stewart [fief. 6 ] made an experiment of the convex 
hull, cheapest angle insertion algorithm (CCA) which will be 
discussed in the next section as a stand-alone algorithm and 
with each of the three post- process ors . The algorithms are 

designated CCA, CCA2, CCAO, and CCA3 for the convex hull 
cheapest insertion stand-alone, with 2-opt, with Oropt and 
with 3-cpt respectively. He drew two conclusions from his 
computational results. First, the 3-opt requires substan- 
tially more time than either the 2-opt or the Oropt. 
Second, the 2-opt is dominated by the Oropt and the 3-opt in 
quality of solution. 

In computation time, Oropt only looks at 

2 

approximately 3n of the n possible 3-opt exchanges on 
each pass. There are n ways to select the first branch, 
times 3 ways to select the second branch, and n-2 ways to 
select the third branch. 

This accounts for the fairly close times for the 
2-opt and Oropt. The quality of CCAO solutions dominate 
CCA2 solutions. On the other hand, there is little or no 
difference between the Oropt and 3-opt in terms of solution 
quality . 

Stewart’s main conclusion from the above experiment 
is that the Oropt performs as well as a 3-opt in a small 
percentage of the computer time required by a 3-opt, and it 
should be preferred to both the 2-opt and the 3-opt for 
Euclidean TSP's. 
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3. Composite Proced ure 



The tasic composite procedure is a combination of 
the tour construction and branch exchange procedures. It is 
obtained by appending a branch exchange procedure to the 
tour construction algorithm as a post- processor. The proce- 
dure can be stated as follows [Ref. 17]. 

Step 1. Obtain an initial tour using one of the 
tour construction procedures. 

Step 2. Apply a branch exchange procedure to the 
solution produced by the step 1. 

Stop wnen no further improvement can be made. 

The composite procedure is relatively fast computa- 
tionally and gives good results [Ref. 18]. 

E. CCAC 



1 - Al g orithm 

The GCAO algorithm designed by Stewart [Ref. 6] uses 
the convex hull of the nodes in N for its initial suntour. 
Then it inserts the nodes not currently in the subtour where 
they may be inserted most cheaply (the Cheapest Insertion 
criterion). It selects the node k to be inserted at each 
iteration according to how large an angle is formed by the 
two arcs that must be added to the current subtour 
(Selection criterion) in order to insert k. Finally it uses 
an Oropt to make local improvement on the tour constructed 
in the first stage. CCAO means Convex Hull, Cheapest 
Insertion, Angle Selection, Oropt. 

Algorithm : CCAO 

Input : Number of nodes, x and y co-ordinates of all 
nodes. 

Output: Ordered list of tour, total cost. 
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step 1 . (Initial Subtour) 

Find the convex hull of the set of nodes N. 
Call the set of nodes on the boundary N'. 
Let the initial subtour be the nodes of N* 
in the same order as they appear on the 
convex hull. 

Step 2 . (Cheapest Insertion) 

For each node m€N-N’, find 

(i , j ) = argmin c + c - c 

m m i,j im mj ij 

s.t. i, j € N*, i, j : connected. 



Step 3 . (Greatest Angle Selection) 

For the next insertion, select the node 

that maximizes the angle between the arcs 

(i , m) and (m,j ) over all m € N-N*. 
m m 

I.e, find k = argmax angle { (i ,m),(m,j ) } 

m ID ID 



s.t. m € N-N’ . 

Insert k between i and 1 

k k 

k to N ’ . 

step 4 : If N* = N, go to step 5- 

Otherwise return to step 2. 
step 5 ; Apply an Oropt to the current 
when no further improvements can 
End of algorithm CCAO 



and add 



tour. Stop 
be found. 



2. Exa mcle 

Figures 2.2 - 2.4 illustrate the above algorithm on 
the TSP defined as test problem £ 1 ] in Appendix A. First 
the convex hull is generated for an initial starting 
subtour. This subtour consists of nodes 2,13,12,14,5,15,7,4- 
A solid line marks the boundary of the convex hull in Figure 
2 - 2 . 
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Figure 2.2 Initial Sabtonr and Issertion. 

In step 2, each of the interior nodes (1,3,6,8,9,10,11,16) 
is associated with a pair of connected nodes on the initial 
subtour (the dashed lines in Figure 2.2). In step 3, the 
dashed lines that form the greatest angle (closest to 18Q0) 
identify the node to be inserted (node 10 in this example). 

Figure 2.3 shows the problem after the first three 
insertions ( node 10, node 1 and node 8 in that order). 
Notice that some nodes not in the subtour are associated 
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Figure 2.3 Internedia te Sabtoar and Insertions. 

with new node pairs. Figure 2.4 shows the final tour for 
stage one. This tour is now passed to an Oropt post- 
processor. In this case the tour from stage one appears 
from inspection to be optimal, and Oropt will find no 
improve irent . 
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X 



C -J 

Figure 2.4 Final Tour of CCIO- 
3. Computational Result s 

In addition to CCAO, CCCO (Convex, Cheapest, 
Cheapest, Oropt) has been coded for the purpose of compar- 
ison. The only difference between CCAO and CCCO is that 
CCCO uses the cheapest selection criterion instead of 
greatest angle of CCAO. 

Me used Sedgewick’s [Ref. 19] package wrapping algo- 
rithm for finding the convex hull (initial subtour). 
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starting with soaie point (called the anchor) that is guaran- 
teed to be on the convex hull (say the one with the smallest 
y co-ordinate), take a horizontal ray in the positive direc- 
tion and sweep it'upward until hitting another point. This 
point is on the hull. Then steirt at that point and continue 
sweeping until hitting another point, etc. The package is 
completely wrapped when the first point is included again. 
The following algorithm finds the convex hull of an array 
l(1,...,n) of nodes, the node L(n+1) is used as a sentinel, 
that is, a copy of the first node which is used to signal 
completion of the procedure. The variable NH is maintained 
as the number of nodes so far included on the hull- 

hlgorithm i Package trapping 

Input : Number of nodes, x and y co-ordinates or all 
nodes. 

Output: Ordered list of convex hull and number of 
nodes included on the convex huil- 

£tep 1 . (Initialization) 

find and duplicate anchor. i.e., find 

NaiN = argmin y s. t. i € N and set 

i 

NH = 0 , 1 (n + 1) = L (NMIN) . 

Step 2 ; (Swap nodes NH and NHIN), 

Put last node found into the hull by 
exchanging it with the NHth node. 

NH = NH ♦ 1. 

TEHP = L (NH) . 

L(NH) = L (NM IN) . 

L(NHIN) = TENP. 

Step 3 : (Compute angle) 

Compute the angle from the horizontal made 
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step 4 



by the line between L(NH) and each of the 
nodes not yet included on the hull. 

: (Find next hull node) 

Find the node whose angle is smallest among 
those with angles bigger than the current 
value of the ‘sweep* angle (the angle from 
horizontal to the line between L(NH-1) 
and L(NH) ) . 

Step 4 : Stop when the first point is encountered 
again. I.e.# L(n+1) = L(NililN) . 

Otherwise, go to step 2. 

End of algorithm Package Wrapping 

He used Sedgewick's Pseudo Angle for finding the 
smallest angle in step 3, which is coded as the 'THETA* 
function. This function returns a real number between 0.0 to 
4.0 that is net the angle made by LI and L2 with the hori- 
zontal but which has the same order properties as the true 
angle. If dx and dy are the delta x and y distances from 
some node to the anchor node, then the angle needed in this 
algorithm is arctangent dy/dx. However, the arctangent 
function is likely to be slow and it leads to at least two 
annoying extra conditions to compute : whether dx is zero, 

and which quadrant the point is in. 

In this algorithm we only need to be able tc compare 
angles, not measure them. Thus it makes sense to use a func- 
tion that is much easier to compute than the true angle but 
has the same erdering properties as the true angle. A good 
candidate for such a function is simply dy / (dy ••• dx) . 
Testing for exceptional conditions are still necessary, but 
simpler. 

Function THETA ( Pseudo Angle ) 

Input : dx,dy (delta x and y distances from some 
node to the anchor node) . 
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Output : Pseudo angle made by II and L2 with the 
horizontal line. 



begin 








dx = 


X (12) 


- X (LI ) : ax = 


abs(ax) : 


dy = 


y (12) 


- y(L1) : ay = 


abs (ay) : 


if ( 


dx=0 ) 


and ( dy=0 ) 


th en t = 0.0 








else t = dy / (ax + 


if dx 


< 0 


then t = 2.0 ■ 


- t 






else if dy < 0 


then t = 4.0 + t 


end 









End o£ function THETA 

Figure 2.2 shows how the hull is discovered in this 
way. We used Sedgewick's pseudo Angle for finding the 
greatest angle selection point. 

The data for our test problems is given in the 
Appenaix. The computational results are summarized in Table 
I. As can be seen in Tacle I, CCAO is faster than CCCO on 
the small-scaled test problems (below 30 nodes ), but CCCO 
is faster than CCAO on the moderately large sized problems 
(over 5 0 nodes) . Generally, the accuracy is almost iden- 
tical in both cases. 

Stewart £Ref. 6] showed that in a large scaled 
problem, the CCAO algorithm outperforms any other insertion 
and selection algorithms. Thus, we are highly motivated to 
use a modification of the CCAO algorithm for solving the 
time-window ccnstrained TSP. 
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TABLE I 

CCBBOTATIOHAL BESniXS OF CCCO, CCAO 



CCCO CCAO 



Prob lem 
Number 


Number 

of 

Nodes n 


Best 

Known 

Solution 


5 ; 

Over 

Best 


CPU 

Time 

(sec) 


cr 

/o 

Over 

Best 


CPU 
. Time 
(sec) 


7 n 


16 


66.603 9 


0.00 


0.0133 


0.00 


0. 00t>6 


[2] 


22 


469. 0288 


0.00 


0.0233 


0.00 


0.0100 


[3] 


22 


278. 4371 


0.00 


0.0166 


0.00 


0.0066 


[5] 


51 


429. 7000 


2.72 


0. 1897 


3. 94 


0.2562 


[ 6 ] 


76 


552.9000 


1. 64 


0.5857 


1.54 


0-6889 



* CPU times it seconds on IBM 3033 



III. I EE ISP WITH HARD TIHE WINDOW CONSTRAINTS 



A. INTBODUCIION 

The first time-constrained TSP we consider is the case 
in which late arrivals are not allowed, and early arrivals 
must wait for the opening of the time window before they can 
begin to service a customer. This is called the hard time 
window case and it is illustrated in Figure 3.1. 



cannot be 
violated 



waiting 

— : — >1 time |< 

X <- 

salesman lower 

arrives bound 




time window 
for city i 



— > 

upper 

bound 

1 ^ 



. Figure 3. 1 Diagram for Hard Time Window Case. 

The hard time window case corresponds to military opera- 
tions and to some civilian distribution problems. Meeting a 
deadline is considered a critical factor in this case. The 
soft time window case will be .discussed in the nextchapter. 

Consider a graph G = {N,A} composed of a set of nodes N 
and a set of arcs A connecting these nodes. We now define 
some notation to be used throughout cur discussion of the 
time-window-ccnstrained TSP. 
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1 = Lower bound on the time window at node i 

(early allowable arrival time at city i) . 

u = Opper bound on the time window at node i 

(latest allowable arrival time at city i) . 

d = Time required to spend at node i. 
i 

(service time at city i) . 

SPEED = Constant speed at which the vehicle 
travels. 

dist = Distance from i to j. 

i j 

c = Travel time from i to j. 

ij 

Note ; c = dist / SPEED. 

ij ij 

We use c and c(i, j) interchangeably. 

ij 

depot = Depot (home) node. 

L = ( L(1) ,L(2) ,...,L(n) ). 

= A tour with n stops visited in the order 

l(^)»l(2)»...»L(n). 

ARBVT = Arrival time at city i. 
i 

WAIT = Waiting time at node i for the hard time 
i 

window. 



He also use l(i)/ u(i), d(i),. ARRVT(i), WAIT(i) and 

1 , u , d , ARBVT , WAIT interchangeably, 
i i i i i 



B. HEORISTIC SOIUTIOM TECHHIQOES FOB BARD TIME HINDOHS 
1. Ne a r est Neighbor 

The following is a Nearest Neighbor heuristic 
similar to the one used in the unconstrained TSP. At each 
iteration we add a new node to the end of the subtouc. It 
is the first node that can be visited from the last node of 
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the subtour, taking into account any waiting time that might 
te necessary due to the lower time window bounds. 



Algorithm : Nearest neighbor 

Input : Number of nodes, x and y co-ordinates of all 
nodes, time windows for all nodes. 

Output : Ordered list of tour, total travel time. 



Step 1 
Step 2 

Step 3 
Step 4 

S tep 5 

Step 6 

Step 7 
Step 8 



(Initialization) 

Start at the depot. 

Let i=depot, N‘ = {i} . 

Compute AERVT for all nodes k€ N-N’ if k can 

te visited directly after i ; 

AERVT = max { AERVT , 1 } + d + c 

k i i i ik 

If ARRVT > u , then stop (• no feasible 
k k 

solution* ) . 

If ARRVT < 1 , then cost = 1 . 

k k k k 

Otherwise, cost = ARRVT 

k k 

( Nearest Neighbor Selection) 

Choose the node k€ N-N* such that cost is 

k 

a minimum. I.e, find 
k = argmin cost s.t. j €N-N*. 

j 

(Insertion) 

Insert k after i, add k to subtour N*, and 
let i = k. 

If N’ = N, go to next step. 

Otnerwise, return to step 2. 

Compute total travel time, then stop. 

Total travel time = max { ARRVT , 1 } 

k k 

d + c 

k k, depot 



End of algorithm Nearest Neighbor 
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Tills solution was constructed by starting at the 
depot and moving to the nearest neighboring customer that 
has not yet teen visited as long as the upper bound level 
was not violated. This heuristic may fail to solve the 
problem . 



2. SCCO 

This algorithm is designed for the case when some of 
the nodes do not have time windows. Se call these nodes •' 
time free 

SCCO is similar to the cheapest selection, cheapest 
insertion method for the unconstrained TSP, except that the 
nodes with time windows are treated differently from the 
time free nodes. The nodes with windows are inserted in 
order of increasing upper time window bound. 

The time free nodes are inserted between these nodes 
by cheapest selection and cheapest insertion, for as long as 
the upper bound time window will allow. in the end, a 
Modified Oropt is used to improve the solution. 

There is one possible difficulty with this approach. 
It may become impossible to reach some of the time- 
constrained nodes before their upper bound. In this case, 
we must delete some node (s) from the subtour. Whenever we 
see an upper nound that cannot be satisfied, we select a 
node to delete by the following criteria. 

The first criterion is the width of the time window. 
Hence, time-free nodes are considered first. Then, if 
several nodes in the subtour are tied for the widest time 
window, we select for deletion the node that results in the 
greatest time saved. The algorithm is summarized as 
follows . 

Algorithm : Successive Cheapest Cheapest Oropt (SCCO) 

Input : number of nodes, x and y co-ordinates of ail 



35 



0 utput 



nodes, time windows for all nodes, 
ordered list of tour, total travel time. 



Step 1 . (Initialization) 

Start at the depot. 

Let i=depot, N* = {i}. 

Set k = argmin u s.t. j € N-N'. 

j 

If k is time free node, then set k = depot. 

Calculate AERVT . 

k 

AHRVT = max { ARRVT , 1 } + d + c 

k i i i ik 

If AERVT < u , then go to step 5. 
k k 

Otherwise, select time free node m € N* which 
results in the greatest time saved for dele- 
tion. Delete node m from N’, go to step 3. 
Step 5 . Add node k to the subtour N'. 

Step 6 . Insert time free node j€ N-N* between nodes i 
and k by cheapest insertion and cheapest 

selection (same as CCCO) until AERVT does not 

k 

exceed u 

k 

If ARRVT < i , then set ARRVT = 1 . 

k k • k k 



Step 2 . 
Step 3 . 

Step 4 . 



Step 7 . If N'=N, then go to next step. 

Otherwise, let i = k and go to step 2. 

Step 8 . Apply the Modified Oropt procedure to the 
current tour. Stop when no further improve- 
ments can be found. 



End of algorithm SCCO 

"Successive" means select the node successively by 
the smallest upper bound. In the SCCO algorithm, if the 
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salesman arrives before tie lower bound of the time window, 
adding waiting time, we set the arrival time equal to the 
lower bound. 

The flcdified oropt procedure for- improving the solu- 
tion is described below. This procedure consider only those 
exchanges that would result in a node being inserted between 
two other nodes in the current tour. 




Figure 3.2 Current Tour before Hodified-Oropt. 

Figures 3-2 and 3.3 are helpful to understand how 
the procedure works. In both figures, i,j,k,l, and m are the 
nodes in the current tour. Nodes 1 and m are considered to 
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be adjacent tc node Jc. A test is then conducted to deter- 
mine if node k can be located between two other nodes, such 
as i and j, sc that it results in reduced total travel time- 
l£ it can, we make the appropriate arc exchanges, then 
update the total cost and route orders. 




In this example, the three arcs (i, j) , 
<k,m) are removed and replaced by (i,k) , (j#k). 

When no further exchanges improve the solution, 
rithm terminates. 



^k ,1) , and 
and ( 1, m) . 
the algo- 
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Figure 3.4 Subtoar in SCCO Procedure. 

Figures 3. 4 - 3.6 illustrate the SCCO algorithm for 
the TSP with hard time windows given in Appendix F as in 
test problem [1]- In this problem 10 of 16 nodes have time 
windows- The ether 6 nodes are time free. 

First, the subtour starts at the depot (node 16) and 
we insert the node with the smallest upper bound (node 12). 
He examine all nodes which could be inserted between 16 and 
12 as long as the upper bound on node 12 is observed. In 
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Figure 3.5 internediata Subtoar in SCCO Procedure. 

this case there is no such node. Then we select the next 
smallest upper bound (node 14), add it to the tour, look for 
nodes to insert before it, and continue in this manner. Now 
we have formed the partial tour 16, 12, 14, 11, 6, 3 as 

shown in Figure 3.4. 

As shewn in Figure 3.5, we can insert 3 time free 
nodes between node 6 and node 3. These insertions are made 
according to the cheapest insertion and cheapest selection 
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Figure 3.6 Final Subtour for SCCO. 

criteria. He do not maJce any further insertions because 
they would cause a time window violation at node 3. 

Figure 3.6 shows the final tour for the SCCO 
heuristic. This tour is passed to a Modified Oropt, but in 
this case it will find no improvement- 
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3. SCAO 



This heuristic is identical tc SCCO except for the 
use of the greatest angle selection criterion for the time- 
free nodes, instead of cheapest selection- 

llgorithm : Successive Cheapest Angle Oropt (SCAO) 



Input : Numher of nodes, x and y co-ordinates of all 
nodes, time windows for all nodes. 

Output : Ordered list of tour, total travel time- 



step 1 



Step 2 
Step 3 



S tep 4 



(Initi aliz ati on) 

Start at the depot. 

Let i=depot, N* = {i}. 

Set k = argmin u s.t. j € N-N'. 

j 

If k is time free node, set k = depot- 

Calculate AflEVT . 

k 

AEEVT = max { ARRVT , 1 } + d + c 

k .i i i ik 

If AivEVT < u , then go to step 5. 
k k 



Step 5 
S tep 6 



Otherwise, select time free node m € N' which 
results in the greatest time saved for dele- 
tion. Delete node m from N*, go to step 3. 

Add node k to the subtour N'- 

Insert time free node j € N-N' between nodes i 
and k by cheapest insertion and greatest 
angle selection (same as CCAO) until ARRVT 

k 



does not exceed u 

k 

If AEBVT < 1 , then set ARRVT = 1 - 
k k k k 



Step 7 .Leti = k. 

If N' = N, then go to next step. 
Otherwise, go to step 2. 
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step 8 . Apply the Modified Oropt procedure to the 
current tour. Stop when no further improve- 
ment can be found. 

End of algorithm SCAO 

This algorithm is same as SCCO except greatest angle 
selection is used instead of cheapest selection, as in SCCO. 

4. SLACK 

This heuristic was suggested by Professor Rosenthal. 
It is designed for the case when the widths between the 
upper and lower bounds of the time windows are relatively 
large. 

In this heuristic, the SLACK is the most important 
concept. The SLACK (i) can be defined as the maximum amount 
of time by which arrival at node i can be delayed 

without causing an upper bound to oe violated for a node 
currently on the tour. 

The SLACK function can be defined as a recursive 
function as follows. 

SLACK(L<i)) = min { u(L(i)) - AERVT(L(i)) , 

SLACK (L(i+ 1) ) +WAIT(L(i)) } 

where 

WAIT(L(i)) = max {0, l(L(i)J - AREVT (L (i) ) } 

The first element of this recursive function is the 
difference between the upper bound and arrival time at node 
L(i). The second one is the sum of next node's SLACK and 
waiting time of node L (i) . The minimum of these two 
elements is a possible delay time of the arrival time at 
node L(i) without violating the upper bound of ail nodes 
after L (i) in the current tour. 
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The advantage of this recursive function is that it 
is easy to calculate a possible delay time without calcu- 
lating new arrival times for all nodes after L (i) - The algo- 
rithm is summarized as follows. 



Algorithm : SLACK 



Input : Number of nodes, x and y co-ordinates of all 
nodes, time windows for all nodes. 

Output ; Ordered list of tour, total travel time. 



Step 1 
S tep 2 

S tep 3 
Step U 

Step 5 



S tep 6 



(Initiaiizati on) 

Start at the depot. Let N* = {depot}. 

Sort the upper time windows. 

U = ( u ,u ,. . .,u ) 

1 2 n 

s.t. u<u<...<u 
1 2 n 

Set k = argmin u s.t. i€ N-N*. 

i 

Find a node L (ISTAR) after waich node k 
can be inserted in the current sequence, 
if such a node exists. Go to step 7. 

(The criteria by which we determine if an 
insertion can be made are given below.) 

If there is no such place to insert node k, 
then try to find a node L(ISWAP) in the 
current sequence such that k can replace 
L(lSWAP) and I(ISWAP) has a good chance of 
being reinserted somewnere else. 

Select ISWAP which has the largest time 
window width among candidates for ISWAP. 

If there is no candidate, then stop. 

( * no feasible solution * ) 

Do swap ( add k to M', and delete L(ISWAP) 
from N*, and set k = L(ISWAP) ), 
then update slack and arrival times. 

Go to step 3. 
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step 7 . Select the node which results in the minimujr 
additional travel, i .6, the node Jc which 
minimizes the following quantity. 
c(L(I),k) +c(k,L(I+D) - c(L(I) ,L(1+1) ) . 

Step 8 . Insert k after L(ISTAS) , and add k to N ’ , 
and update slack and arrival times. 

Step 9 . If N*=N, then stop. 

Otherwise, go to step 3. 

End of algorithm SLACK 

Ihis procedure starts with sorting an array 

u , u ,...,u into ascending order using a heapsort 

1 2 n 

[Ref. 20]. Ihis u array is used to select a node k in 
ascending order for insertion. Since the upper bound cannot 
be violated, this step is performed. Then find a node 1(1) 
after which node k can be inserted in the current sequence, 
if a such a node exists. 

There are two tests which must be administered to 
determine if k can be inserted after L(I). First, the 
arrival time at node k if k succeeds L(I), whicn is called 
TZST1 must not he greater than the upper bound u . Second, 
if k precedes L (1+ 1) , then the resulting delay in arrival at 
1(1+1), which is called TEST2 , must not greater than 
SLACK (1 (1+1 )) . We can calculate TESI1, T£ST2 as follows. 

TEST1 = Arrival time at node k if k succeed 1 (i) . 

= max {AERVI <1 (i) ) , l(l(i))} + d(L(i)) + c(L(i),k). 

TEST2 = Delay in arrival at L(i + 1) if k precedes L(i+1). 

= max { TEST1, 1 (k) } + d (k ) ’ + c(k,L(i+1)). 
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If there exists more than one node L(I) after which 
node k can be inserted, we select L (I) according to the 
criterion of least additional travel time. This additional 
travel time, called TEST3, is given by 

TEST3 = c(L(I),k) + C(k,L(I + 1)) + C(1(I) ,L(I+1) ) . 

When we insert node k after L (I) , we update the 
arrival times and SLACKS. In the updating process, we 
compute updated values of SLACK only for the nodes whose 
SLACK actually changes as a result of the insertion. 

If there is no place to insert node k, we call a 
subroutine called 'TSWAP'. TSKAP tries to find a node 
L(ISWAP) in the current seguence such that k can replace 
L(ISwAP) and l(ISWAP) has a good chance of being reinserted 
somewhere else. TSWAP uses TESTl ans IEST2 to find a candi- 
date for ISWAf and then uses a largest time window width to 
select ISWAP. If there exists such a ISWAP, then we do the 
swap and update SLACKS and arrival times and try to insert 
again. 

C. EXACT SOXbTION TECHIIQOES FOB HABD TIME WINDOWS 
1. S tate-Space R elaxati on Procedu re 

A dynamic programming model of the time-constrained 
TSP has been developed by Christofides et al. [Ref. 5]. We 
applied their approach to compute bounding information 
within a branch and bound algorithm. 

Consider the TSP defined on the graph G = {N,A} with 
the time window constraints, where N is a set of all nodes 

of G, and A is a set of arcs. Let fi(j) be the set of nodes 
from which it is possible to go directly to node j. We can 
initially set B(j) = N - { i | Ij + dj^ + c^- > uj } , because 

it is impossible to go directly from ncde i to node jif the 
earlist possible arrival time at node j exceeds the upper 
time window cf node j. 
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Let f(S,j) be the duration of the least time path 
starting at node 1 passing through every node of SS S*=N-{1} 
and finishing at node j. For a given S and j, we can calcu- 
late a minimun arrival time in node j as 



T(S,j) =. min [ f(S-j,i) + d.+ c..]. 

re(S-gflfia)) i xj 



Then, 



f(S,j) = T(S,j), 


if 


1 . < T (S 

“j 


,j) ^ 


= 1 . 

3 

= 00 , 


if 


T (S, j) 


< 1 


if 


T (S, j) 


> 



with the initialization; 



= C (1, j) 


# if 


1 < 
j 


c < u 


= 1 


» if 


c < 


1 


j 


» if 


1j ^ 


j 


— 00 


c > 

1j 


u 

j 



(3.1) 



In equation (3.1) the minimum arrival time in node j passing 
through the nodes in the set S can be described as the sum 
of three terms ; the first is the duration of the least time 
path passing through the nodes in the set S- {j} and ending 
in node i, the second is the time required to spend in node 
i, and the third is the travel time from node i to node j. 

The f(S,j) can be calculated for all subsets S of S’ 
and for all nodes j by using equation (3.1) recursively. 
Finally, the optimum solution can be calculated as 



min [ f(S',i) + d + c ]. 

i e S' i i1 



Since the computer storage requirements increase 
exponentially with the size of the problem, this methcd is 
limited to small problems. jT^e total number of f(S,j), when 
S contains k nodes, is k( j, since f (S,j) must be calcu- 
lated for all subsets S or S', and since each node in S must 
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be considered as a possible end-node j. Therefore the 
storage requirement for f (S, j) in a n node problem, is given 
by [ Ref . 21]. 




(n-1) 



n-2 

2 . 



(3,2) 



The storage requirements to solve a 22 node problem 
exceed 22,020,096- For relaxing this limitation, 

Christofides et al. £Eef. 5] proposed a state space relax- 
ation procedure which is analogous to Lagrangean relaxation 
[Ref. 22] in integer programming. The state space associ- 
ated with a given dynamic programming recursion is relaxed 
in such a way that the solution to the relaxed recursion 
provides a bound which could be embedded in a general branch 
and bound method [Ref. 23]. We describe Christofides et al's 
method for doing this below. 

Consider the dynamic programming formulation (3.1) 
The state variable in that formulation is (S, j) , and the 
stage is the cardinality of S. Let g (S) be a mapping from 
the domain of (S,j) to some other vector space (g(S),j)^ 
Let : 



H(g(S),j) = { (g(S-j),i) | i. e (S-j n r ( j) ) } (3.3) 

Since we are interested in lower bounds to the TSP 

with time constraints, H(g(S),j) in (3-3) may be replaced by 

any larger set that is easier to compute. Thus, H(g(S),j) 

can be defined by the following equation: 



H(g(S),j) = { (g(S-j),i) I i£E(g(S),j) } (3-4) 

where (S- j 0 E ( j) ) i E (g (S ) , j ) . 

Fcr calculating the lower bound of, the problem, equation 
(3.1) can be changed to the following equation; 



T <g(S) , j) .min 

(g(S-g) , 1 ) 



ff (g(S- j) ,i) 
H (g (sf ,]) 



c. .] 

13 



(3.5) 
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= min [f(g(S-j) ,i) + d + c ] (3-6) 

i£E(g(S),j) i ij 



This gives us: 



f(g(S),j) = T(g(S) ,j) , if 1 .< T(g(S), j) <u. 

, if T?g(S) ,j) < 1 . ^ 



= 1 



With the initialization: 



if T (g(S) , j) > u 

3 



f(g(j)»j)= c , if l.<c <u. 

= 1 , xf c <1 

3 .. 13^3 

= 00 , xf C > U 

Ij 3 

Finally, tne optimum solution can be calculated as 

min [ f(g(S’)ri) + d ♦ c ]. 

ieE(g(li),1) i i1 

The mapping can be selected from any separable func- 
tion. Christofides et al. used the following mapping 
function. 

g(S) = jSj. (3.7) 



Then eguation (3.6) becomes ; 

T(lc,j) = min [ f( k-1,i ) ♦ d + c ] 

i e E (k, j) i i j 

where k = | S | > 1 . 

This gives us; 

f(k,j) = I(k,j), if 1 < T(k,j) < u 

= 1 . , if T^k,j) < 1 . ^ 

= 00^ , if T (k, j) > u^. 

3 



(3.8) 



With the initialization; 



f(1,3) = 



c(1, j) 
1 



— 00 



if 1 < c < u 

i( c' < ’5 ’ 

if 

'i j 
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Finally, the optimum solution can be calculated as 

min [ f(JS*|,i) + d + c ]. 

i e E ( 1 N I, 1) i i1 



2. Addit ional Condition 

In the previous section, we discussed Christofides 
et al. * s state space relaxation procedure which provides a 
lower bound on the TSP by reducing a state space in dynamic 
programming. This lower bound is effective in branch and 
bound only if it is a tight bound. This is similar to the 
case in integer programming where the effectiveness of 
Lagrangeaa relaxation in producing bounds is relative to the 
integer programming formulation. A redundant state-space 
condition can be helpful to get a better bound- For this 
purpose, an additional condition was used by Christofides et 
al. to avoid loops formed by three consecutive nodes [Ref. 
5]. This can be done in the following way. 

Let k = ISl. Let f (k,j,1) be the duration of the 
least time path from tne initial state to state (k,j) 
without loops formed by three consecutive nodes. Let 
f(k,j,2) be the duration of the second least time path from 
the initial state to state (k,j) without loops formed by 
three consecutive nodes. Let p(k,j,m) be the predecessor of 
j on the path corresponding to f(k,j,m). With the above 
definition, recursion (3.8) becomes; 

T(k,j, 1) = min [ f(k-1,i,m) + d + c ], (3.9) 

i e E(k, jj i ij 

where m = 1, if p(k-1,i, 1)^j 

=2, otherwise. 

This gives us: 

f(k,j,1) = I(k,j,1), if 1. < T(k,j,1) < u. (3.10) 

= 1. ^ if tU,J,1) < 1. ^ 

3 3 
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= 00 » if T(k,j,1) > u . 

j 

Hecursion for f(k, j,2) car be written in the 
fcllowicg way. 



Let : 



T(k,j,2) = min [ f(k-1,i,m) + d + c ], 



i e E ( Jc , 

i 

where m = 1, if p(k-1,i, 1)/j 
=2, ctherwise. 

This gives us: 

f(k,j,2) = I(k,j,2), if 1. < T<k,j,2) < u. 

= 1. , if T]k,j,2) <1. ^ 

= oP , if T(k,j,2) > u^. 

j 

The initialization is 

f(1,i,l) = c(1,i) , if 1 < c < u 

i 1i i 

= 1 , if c <1 

i li i 

= 00 - if c > u 

li i 



and 

f(1,i,2) = =» 

Finally, the optimum solution can be calculated as 

. min [ f(|S'lri»1) + d + c ]. 

1 £ E(j N| ,1 ) 1 i1 



(3.11) 



(3.12) 



(3,13) 



(3.14) 

(3.15) 



Since the additional condition can avoid consxdera- 
tion of a useful lower bound, we considered f(k-1,i,2) in 
recursion (3.9) and (3.11) only when the predecessor of i on 
the path corresponding to f(k-1,i,1) is j. If we do not 
consider the second least time path in case of ? (k- 1, i, 1) = j, 
then f(g(S),j) dees not guarantee the lower bound of f(S,j). 

For this example, let's consider a 4 node TSP with 
time constraints. Node A is the starting node. D is the 
time free node. The lower bound of node B is 9, the upper 
bound of node B is 11, the lower bound of node C is 19, and 
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the upper board of node C is 21. Suppose service time at 
each node is zero. Figure 3.7 shows an optimal route for 
this problem. 




From equation (3.13) we can get; 
f (1,B, 1) = 10, 

f (1,C,1) = 19, 

f (1,D, 1) = 7.07 



Now applying equation (3.9) recursively with i=1, for k = 2 we 
can get; 
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f (2,B,1) 



00 r 



f(2,C,1) = 19, 

£ (2,D,1) = 17.07 

Similarly, for k=3 

f (3,B, 1) = c« , 
f (3,C,1) = ~ , 
f (3, D, 1) — - 



We can see easily that f(3,D,1) is not a lower bound of 
f ({B,C,C} ,D) . 

3. Bra nch and Bou nd Proced ure 

In this section we introduce branch and bound 
enumeration which is used to eliminate subtours in the solu- 
tion of the state space relaxation procedure. Since the 
state space relaxation procedure is a relaxation of the TSP 
with time constraints, the solution to the state space 
relaxation procedure provides a lower bound on the optimal 
value of the ISP with time constraints. Any heuristic solu- 
tion can provide an upper bound. He denote some notation to 
explain this algorithm as follows. 

FL3D = The lower bound, which is tne optimal solut- 
ion to the state space relaxation procedure, 
on the optimal solution to the TSP with 
time constraints given restrictions at the 
current node. 

Z = Current upper bound. 

STACK = Array which represent decision tree. It con- 
tains arc lists which have the same head in 
optimal tour to the state space relaxation 
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procedure given restrictions at the current 
node. 

£c' ] = Travel time matrix given restrictions at 

ij 

the current node. 

There are two types cf tree search. One is depth- 
first search, the other is breadth-first search [Ref. 24]. 
We used depth-first search since breadth- first search 
reguired substantially more storage. Depth-first search 
simply means that when a separation is defined, one of the 
nodes created by the separation is immediately selected to 
be the next sutproblem, and when a ncde is fathomed, the 
enumeration always backtracks to the most recently created 
live node. 

One of the most important requirements of any branch 
and bound algorithm is tight bounds. The closer the bounds 
are to the optimal solution, the fewer nodes must be enumer- 
ated. We used the SCCO heuristic, which was described in 
section B.2 , as an initial upper bound. The lower bound is 
obtained from equation (3.15). 

To save computing time we need a criterion to decide 
whether or not the branching should be continued. If FLBD 
is greater than Z, then the node is fathomed since explicit 
enumeration need not be extended belcw the current node. 
For branching we consider the arcs which have the same head 

node in the directed graph since each arc must have a 

different head in the TSP solution. If there is' no such 
arc, then that solution is a feasible solution. After all 
nodes of the tree are fathomed, a feasible solution which 

has the same value as the upper bound is an exact solution 

to the TSP with tame window constraints. 

The fcllowing branch and bound algorithm is used in 
the programs written for exact solution. 
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Algorithm : Eramch and Bound Procedure 



Input ; lotal travel time of heuristic, travel time- 
Cutput: Ordered list of tour, total travel time. 



Step 1 

Step 2 
Step 3 

Step 4 

Step 5 
Step 6 



Step 7 



(Initialization) 

Let Z = the optimal solution of SCCO. 
STACK = empty. 



Compute FLBD given restrictions defined by 
[ c' J. If FLED > Z, go to step 5. 



13 

(Construct the tree) 

Put all arc(i,j) which have the same head j 
in directed graph on STACK. 

If there is no such arc, save feasible route 
and update Z = FLBD then go to step 5. 

Let travel time of arc(i,j) which is in the 
top of STACK be infinite, then go to step 2. 
(i.e., c' ==>°.) 

(Backtrack) 

If STACK = empty, go to step 7. 

If travel time of arc(i,j) which is in the 

top of STACK is finite, let travel time of 

that arc(i, j) be infinite, then go to 

step 2. (i.e., c' = °° .) 

i j 

Otherwise, let travel time of that arc(i,j) 
be original travel time of that arc(i,j) and 
remove tnat arc(i,j) from top of the STACK- 
Go to step 5. (i.e., c* = c ) 

ij ij 

(termina txon ) 

If there is a feasible route, then the 
optimal travel time = Z. 

Otherwise, there is no feasible solution. 



End of algorithm Branch and Bound Procedure 
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We present the results of our computational experi- 
ence with the algorithms of this Chapter in Chapter V- 
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IV. THE TSP WITH S OFT IIAE HIND CH CONS TRAINTS 



A. INTRODUCTION 

Tie second time-constrainted TSP we consider is the case 
in which both late and early arrivals are allowed by paying 
a penalty cost. The penalties are allowed to be different 
for early and late arrivals. Tne penalty cost is calculated 
as follows. 

Upper penalty cost = max [ 0/ upper penalty constant 

X ( arrival time - upper bound ) ]. 

Lower penalty cost = max [ 0, lower penalty constant 

X ( lower bound - arrival time ) ]. 

In fact, the upper penalty constant is greater than the 
lower penalty constant in most cases. Figure 4.1 may be 
helpful to understand this case. 

This approach makes every problem feasible, no matter 
what the time windows are, i.e, even if it is infeasible in 
the hard time window case. This reflects a practical ppint 
of view, especially when it is possible to save a great deal 
of mileage by allowing a small amount of time window 
violation. 

In this Chapter, we considered one unit of cost to be 
the same as one unit of time. In real world problems, it is 
possible to get a cost by multiplying traveling time by some 
constant. 

We use the notation Ip and up for the lower and 

k’ k 

upper penalty cost at node k. 
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1 





Figure 1 Diagram for Soft Tiae Hindow Case. 



B. HEOfilSTIC SOiaTIOa TECHIIQOES FOB SOFT TIME iINDOiS 



1- Ne arest Neighbor 



This heuristic is similar to the hard tiae windows 
except it takes into account any penalty cost that might be 
necessary. 

Algorithm : Nearest Neighbor 



I nput 
Output 



Number of nodes, x and y co-ordinates of all 
nodes, time windows for all nodes. 

Ordered list of tour, tctal cost. 



Step 1 . (Initialization) 

Start at the depot. 

Let i=depot, N’ = (i) , cost = 0- 

i 

Step 2 . Compute ARRVT for all-ncdes k € N-N' 

AEEVT = AfifiVr + d ♦ c . 
k i i ik 

cost = cost * d + c 
k i i ik 
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9 



Step 3 



S tep 4 



If ARRVT < 1 
k Jc 

If AERVT > u 
k k 



then cost = cost + Ip 
k k k 

then cost = cost up 
k k k 



( Nearest Neighbor Selection ) 
Select the node k € N-N* such that 
is a minimum. I.e. , find 
k = argmin cost s.t. j € N-N’. 

j 



cost 

k 



Step 5 . ( Insertion ) 

Insert k after i, add k to subtour N', and 
let i = k . 

Step 6 . If N' = N , then go to next step. 

Otherwise, go to step 2. 

Step 7 . Compute total cost, then stop. 

Total cost = cost + d + c 

k k k, depot 



End of algorithm Nearest Neighbor 

This solution was constructed by starting at the 
depot and moving to the nearest neighboring customer that 
has not yet been visited. The term "nearest" is modified in 
the sense that we add a penalty cost to the travel time if 
the time window for city i is violated. 



2. SCCO 

This algorithm is also designed for the case when 
there is a combination of tight time window nodes and time 
free nodes. The strict observance of the upper bound in the 
hard time windows is replaced by a penalty cost. 

Algorithm ; SCCO 

Input : Number of nodes, x and y co-ordinates of ail 
nodes, time windows for all nodes. 

Output : Ordered list of tour, total cost. 

Step 1 , llnitialization) 
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0 . 



step 2. 
Step 3 . 

Step 4 . 
Step 5 . 



Start at the depot. 

let i=depot, N* = {1} , cost = 

i 

Set k = argmin u s.t. j€N-N*. 

If k is time free node, then set k = depot. 

Insert node k in the subtour N'. 

Compute ARRVT 

k 

ARRVT = ARRVT ♦ d + C . 
k i i ik 

Insert time free node j € N-N' between nodes 

i and k by cheapest insertion and cheapest 

selection ( same as CCCC) until ARRVT does 

k 

not exceed u . 

k 

Update cost 

k 

cost = cost + d + c 
k i i ik 

If ARRVT < 1 , then cost = cost + Ip 
k k k k k 

If ARRVT > u , then cost = cost * up 
k k k k 'k 



Step 6 „ let i = k. 

If N* = N, then go to next step. 

Otherwise, go to step 2. 

Step 7 . Apply the Modified Oropt procedure to the 

current tour. Stop when no further improve- 
ments can be found. 



End of algorithm SCCO 

This procedure is also similar to the cheapest 
selection, cheapest insertion method for the unconstrained 
TSP, except that the nodes with time windows are treated 
differently from the time free nodes. The nodes with time 
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windows are inserted in order of increasing upper time 
window bounds. The time free nodes are inserted between 
those nodes by cheapest selection and cheapest insertion, 
for as long as the upper bound of the time windows will 
allow. 

In the end, a Modified oropt is used to improve the 
solution. This procedure consider only those exchanges that 
would result in a node being inserted netween two other 
nodes in the current tour. 



3, SCAO 

This algorithm is also designed for the time window 
set which is composed of some tight time windows and some 
time free nodes- 



Algorithm : SCAO 

Input : Number of nodes, x and y co-ordinates of all 
nodes, time windows for all nodes. 

Output : Ordered list of tour, total cost. 



S tep 1 . 



Step 2. 



Step 3 . 



Step 4 . 



(Initialization) 

Start at the depot. 

let i=depot, M' = (i), cost = 0. 

i 

Set k = argmin u s.t. j€N-N'. 

3 

If k is time free node, set k = depot. 

Insert node k in the subtour N'. 

Compute AfiRVT 

k 

AERVT = AEfiVT + d + c 

k i i ik 

Insert time free node j € N-N* between nodes 

i and k by cheapest insertion ana greatest 

angle (same as CCAO) until AfiEVT does not 

k 

exceed u . 

k 
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step 5 . Update cost 

k 

cost = cost * d ♦ c 
k i i ik 

If ARRVT <1 , then cost = cost + Ip 

k k k k k 

If ARRVT > u , then cost = cost + up . 
k k k k k 

Step 6 , let i = k. 

If N' = N, then go to next stop. 

Otherwise, go to step 2. 

Step 7 . Apply the aodified-Orop t procedure to the 

current tour. Stop when no further ioprove- 
ments can be fcund. 

End of algor itho SCAO 

This algorithm is same as SCCO except a greatest 
angle selection in stead of a cheapest selection in SCCO. 

C. EXACT SOIOTICM lECHHIQOES FOB SOFT TIME iINDOiS 
1. St at e-S pace Relaxati cn Pr o ce dure 

In this section we describe a state space relaxation 
procedure, which is adapted from Christofides et ai. [Ref. 
5], for soft time windows. They only considered the TSP with 
hard time windows and without time windows. The difrerences 
are as follows. Tne waiting cost is replaced by a penalty 
cost to be paid in the early arrival case. Late arrival is 
allowed, but a penalty cost has to be paid. So we have to 
calculate tne duration and the penalty cost on each possible 
path to decide the least cost path in each stage. Ke denote 
the penalty ccst on each possible path as PC in this 
sect ion . 

Consider the TSP defined on the graph G = {N,A} with 
scft time window constraints. Let S* be a set of all nodes 
except starting node. Let S be a subset of S'. Let f (S,j) 



62 



b€ the cost of the least cost path starting at node 1 
passing through every node of S and finishing at node j. 
Let T(S,j) be the total duration of a path corresponding to 
f(S,j). Let p(S,j) be the predecessor of j on the path 
corresponding tc f {S, j) . Let lp(t) be the early arrival 
penalty cost function and up (t) be the late arrival penalty 
cost function. For a given S and j, total duration of a 
path can be calculated as 

T(S,j) = [ T(S-j,i) d. ♦ c. . ]. (4.1) 

where p (S, j) = i. 

In equation (4.1) total duration of the least cost path 
passing through the nodes in the set S and ending in node j 
can be described as the sum of three terms: the first is 

total duration of the least cost path passing through the 
nodes in the set S-{j} and ending in node i, the second is 
the time required to spend in node i, and the third is the 
travel time from node i to node j. The dynamic programming 
recursion to determine the least cost path may then be 
stated as 



f(S,j) = min [f(S-j,i)^-d + c ^PCj 
1 e S- j i i j 



where T1 = 


[ 1(S 


-jri) ♦ . 


+ 




c 


]. 






1 






ij 


PC = 


0 


, if 1. < 


T 


1 


< 


u 


= 


Ip (1 . 

-1 


-T1) , If 


T 


1 


< 


1^ 


= 


up (T 1 


-u ) , if 


T 


1 


> 


3 

u ^ 



with the initialization: 



t ({1} / 1) = C 


, if 1 < 


c < u 




ll 


3 


^3. 3 




= c 


♦ Ip (1 -c ) 


, if c 


< 1 




1 ll^ 


13 


3 


= . c 


+ up (c -u ) 


, If c 


> u 


13 


LJ. 

LJ. 


1j 


3 


Finally, the optimum 


solution can 


be calculated 


min [ f (S* , i) 
i e S' 


> d + c 

i i1 


]. 





(4.2) 
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since the computer storage requirements are 
increased exponentially with the size of the problem, this 
method is limited to small problems. For relaxing this 
limitation, a state space relaxation procedure can be used 
same as Chapter III. 

Consider the dynamic programming formulation <4. 2) 
The state variable in that formulation is (S,j), and the 
stage is the cardinality of S. Let g (S) be a mapping from 

the domain of (S,j) to ■ some other vector space (g(S),j). 

Let : 

H(g(S),j) = { (g(S-j),i) I i€S-j} (4.3) 

Since we are interested in lower bounds to the TSP 
with time constraints, H(g(S),j) in (4.3) may be replaced by 
any larger set that is easier to compute. Thus, H(g(S),j) 
can be defined by the following equation: 

H(g(S),j) = { (g(S-j),i) 1 i€E(g(S),j) } (4.4) 

where £- j £ E (g (S) , j) . 

For calculating the lower bound of the problem, recursion 
(4. 1) can be changed to the following equation: 

I (g(S) , j) = [ T (g(S-j) ,i) * d. * c. . ] (4.5) 

1 ID 

where p(g( S) ,d) = i. 

Recursion (4.2) may be stated as 



f(g(S),j) =. min [fig(S-j> 

(g(S-D) ,1) e H(g(5f,j) 


,i) + d + c 


PC] 


(4-6) 


i ij 




= min Cf (g (S-j) 


, i) + d ♦ c 


* PC] 


(4.7) 


X e E (g (S) , j ) 


i ij 





where T1 = 


L KgCS-j) ,i) d . ♦ 


c. . 1# 




X 


ID 


PC = 


0 , if 1 , < T 1 

Ip a .-T1) , If T 1 

-j 


< U ' 


= 


< 1^ 
-i 


= 


up (T T-u ) , if T 1 


> u^ 




j 


D 
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with th€ initialization: 

f (g (j) » j) = c , if 1 . < c < u . 

? -P 1 

= c * lp(l -c J , If c < 1 . 

Id d Id ^ ^ 

=c +up(c -u),if c >u 

Id ij d Id d 



Finally, the optimum solution can be calculated as 

min rf(g(S’)#i)+<i ]. 

i E E (g (K) , 1> i i 1 

The mapping can be selected from any separable func- 
tion. We used a mapping function (3.7) , which is proposed by 
Christofides et ai. , same as Chapter III. Then equation 
(4.5) becomes: 



T (|S| , j) = [ T(|S|- 1,i) + d. + c. . ]. 

1 13 

where p(|S | , j) = i 

Eecursicn (4.7) may be stated as: 



(4.8) 



, 3 ) =. min [f((S|- 
i e E(lSt ,3) 


Iri) * 


[. 1(1S(-1, 


i) + 


d + 


C. . ] 






1 


ID 


0 » if 


1 . 


< T 1 


< u 


ip (1 -Tl) 


j . 
» ir 


T 1 


< 1^ 


up (T?-u ) 


» if 


T 1 


> u' 



c + PC ] 

ij 



PC = 0 



with the initialization: 



f (1»D) = c 



= c 



Ij 

Ij 



if 



1 < c < u 

j ^ j j 

+lp(l-c ),lf C<1 
= c ♦ up (c -u ) , if c > u 

Ij Ij j Ij j 



(4.9) 



Finally, the optimum solution can be calculated as 



min [ f(lS'i,i) + d + .c ]. 

i e E(|N|,1) i i1 
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2 . 



Additional Condi tion 



In th€ previous section, we discussed a state space 
relaxation procedure which is adapted from Christofides et 
al.fRef. 5]. That procedure provides a lower bound on the 
ISP with soft tine window constraints. The additional 
condition to avoid loops formed by three consecutive nodes 
was used to get a better bound [ Eef , 5]. This can be done 

in the following way. 

Let k = jS|. Let f(k,j,1) be the cost of tne least 
cost path from the initial state to state (k, j) without 
loops formed by three consecutive nodes. Let f(k,j,2) be 
the cost of the second least cost path from the initial 
state to state j) without loops formed oy three consecu- 
tive nodes. Let P(k,j,m) be the predecessor of j on the 
path corresponding to f(k,j,m). With the above definition, 
equation (4.8) becomes; 

T(jc,j,m') = [ T(k-1,i,m) + d + c ], m' = 1,2 (4.10) 

i ij 



where F (k, j , n' ) = i 



m = 1 , if p (k- 1, i, 1 ) # j 
= 2 , otherwise. 

With the initialization; 

T (1, j, 1) = c_ 

1] 



and 



T(1,j,2) = 00 . 

Recursion for f (k,j,l) can be calculated in the following 
way. Let: 

T* (k,j,m) = £T(Jc-1,i,m) + d ♦ c ], m = 1,2 

i ij 

Ibis gives us; 
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f(k, j, 1 ) = min [£(k-l,i,m) + d, + c + PC] (^- 11 ) 

i e E(k, jj 1 xj 

where PC = 0 , if 1 . < T* (k,j,m) < u. 

= Ip (1 -T*(k, j^m)) , if T*(k, j,m? < 1 . 

= up (T^ (k, j,m)-u ), if TMk, j,m) > u. 

j D 

IB = 1 , if P(k-1,i,1)#j 

= 2 , otherwise. 



. With the initialization: 



f (1,i, 1) = c^. 

1i 

= c 

li 

= c 

li 



, if < 

+ lp(l.-c^.) 

X lx 
+ up(c -u ) 
1 i i 



c . < u . 



li 


i 


if 


c <1 




li i 


if 


c > u 




li i 



(4.12) 



Recursion for f (k, j,2) can be written in the following way: 



f (K, j , 2 ) 


= min rf(k-1,i,m) * d + 


c 


ieE(k,J) X 

X5tp(k,g, 1) 


ij 


where PC = 0 


, if 1, < T» (k, j,m) < u . 
-T* (k, j^m)) , if T* (k, j,m? 




= lp(l 


< 1 


= up (T 


^ (k, j ,m)-u .) , if T’(k, j,m) 


> u 



— — 'S'- W ^ W — W — ^ W ^ w w , 

1 1 

m = 1 , if p(k-1,i,1) 

= 2 , otherwise. 



(4. 13) 



With the initialization: 

f(1,i,2) = oo (4.14) 

Finally, the optimum solution can be calculated as 



min [f(|S'|/ir1) 
e E(|N1,1) ' 



+ c 



i1 



( 4 . 15 ) 



Since the additional condition can a7oid considera- 
tion of a useful lower bound, we considered f(k-1,i,2) in 
recursion (4.11) and (4.13) only when the predecessor of i 
on the path ccrres pondinq to f(k-1,i,1) is j. If we do not 
consider the second least cost path in case of p (k-1 ,i, 1) = j , 
then f(g(S) , j) does not guarantee the lower bound of f(S,j) 
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For this example, let's consider 4 node TSP with 
time constraints. Node A is the starting node- D is the 

time free node. The lower bound of node B is 9, the upper 
bound of node B is 11, the lower bound of node C is 19, and 
the upper bound of node C is 21. Suppose service time at 
each node is zero, lp(t)=t, and up(t) = 5t. Figure 3.7 shows 
an optimal route for this problem. Frcm equation (4.10) and 
(4. 12) we can get: 



f(1,B,l) = 10, 1(1, B,l) = 10, p(1,B,1) = A; 

f(l,C,l) = 19, T(1,C,1) = 14.14, p(l,C,l) = A; 

f(1,D,1) = 7.07, T(1,D,1) = 7.07, p(1,D,1) = A. 

New applying equation (4.10) and (4.11) recursively with 
i=1, for k= 2 we can get: 

f (2,3,1) = min [ 94.7, 29.84 ] = 29.84, 
i s (C,D) 

T(2,3,1) = 14.14, p(2,B,1) = D. 

Simila rly, 

f(2,C,1) = 19, T(2,C,1) = 14.14, p(2,C,l) = D; 

f(2,D,l) = 17.07, T(2,D,1) = 17.07, p(2,D,l) = B. 

F cr k = 3 , 



f (3,B, 1) 
T (3,B, 1) 



= min [ 
i e{C} 

= 24. 14 



94.7 ] = 94. 7, 

, p(3,B,1) = C. 



Similarly, 



f(3,C,l) = 39.84, T(3,C,1) = 24.14, p(3,C,l) = D; 
f (3,D, 1 ) = ^ . 



We can see easily that f(3,D,l) is not a lower bound of 
f ({B,C,D} ,D) . 
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3. Branc h and Bound frocedure 



We used the same branch and bound procedure used to 
eliminate subtours in the solution of the state space relax- 
ation procedure in Chapter III-C. 3. The SCCO heuristic, 
which was described in section B-2, was used as an initial 
upper bcund, and the lower bound was obtained from equation 
(il. 15) . 

We present the results of our computational experi- 
ence with the algorithms of this Chapter in the next 
Chapter . 
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V. COMPDTATIOKAL EIP EBIEN CE 



1. TESI PROBIEHS 

Four sets of test data are used in this thesis. lest 
problem number [1] is taken from Sedgewick £Ref.19: p. 309]. 

The other problems, numbered [2], [3] and [4], are from 

Appendix 9,1 of Eilon et al. • s text [Ref. 21]- These test 
problems are shown in Appendices A,B,C respectively. These 
published problems contain node and depot locations, but 
they do not include time windows. 

We constructed time windows for test problems 
[1]»f2],[3] by first using the CCAO neuristic on the uncon- 
strained TSP. lime windows were then placed aoout each node 
such that the CCAO route was feasible. The idea for gener- 
ating time windows in this way comes from Baker [Ref. 25], 
who used the unconstrained Nearest Neighbor heuristic as his 
starting point instead of CCAO. 

The time window widths were set to varying sizes ranging 
from 3 to 14. Some of the time windows were fairly tight 
while others overlapped. This is in contrast to Baker’s 
work, where all the time windows have width equal 2 units. 

The last problem number [4] is the same as test problem 
[3], except that the time windows were constructed from a 
Nearest Neighbor solution to the unconstrained traveling 
salesman problem, as in Baker [Ref. 25]. Figure 5.1 

displays the CCAO solution for test problem [3] and Figure 
5.2 illustrates the unconstrained Nearest Neighbor solution 
for the test problem £4]. We found a small error in Baker’s 
TSP solution for the Nearest Neighbor [Ret. 25], in that the 
nearest node from node 16 is node 17, not node 13. The 
resulting cost is actually higher, it is 312.09, not 310-22. 
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Figure 5.1 Unconstrained Solution Obtained by CCAO. 

Each of the four sets of test data was used to create 
four test problems. The separate instances differed in the 
percentages of time window constraints that were chosen to 
be in effect. The four cases were 100%, 90%, 75%, and 50%. 
We refer to this percentage as the "time-window percentage". 

A random number generator was used to decide which nodes 
would have time windows. Test problems for the time window 
constrained TSP are shown in Appendices G through V. The 
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Figure 5.2 Dncon strained Solution Obtained by 
Nearest Neighbor Heuristic. 

penalty cost factors can be varied depending on real world 
problems. We used 2 and 5 as the lower and upper penalty 
cost factor. Also we set th a service time at each node to 0 
to make it easy to construct the time windows. 

The computational results are presented in Tables II and 
III. The figures reported represent results of our test 
runs for each test case. 
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B. COBfOTATIONAX BESOLIS 



I* Hard lime Windows 

As noted in section 3, the Nearest Neighbor 
heuristic often cannot solve the problem, because the 
arrival time of the nearest node frequently violates the 
upper bound. However in test problem [4], the results of 
the Nearest Neighbor are the same as in the unconstrained 
problem, because this problem itself was constructed by a 
Nearest Neighbor heuristic. 

Generally, the SCCO and SCAO heuristic can be easily 
applied to the hard time window TSP. According to our 
experiments, if the time window width becomes large relative 
to the travel time between nodes on the optimal uncon- 
strained TSP route, then the lower percentage time window 
problems become more difficult to satisfy. This phenomenon 
can be seen in test problem [1]. That is because the other 
nodes in the optimal route for the unconstrained TSP problem 
could be inserted without causing violation of the upper 
bound. 

The SIACK heuristic takes slightly more time than 
the other heuristics. It achieved lower accuracy in test 
proolems [1 ] and [3] in the 50 percent time window case. 
The exact algorithm can find the exact answer in most prob- 
lems, but when there are fewer windows in effect, it takes 
more computation time. It cannot solve the 50 percent time 
window problems [2] and [4] within 180 seconds. 



73 



COHPOTATIONAL RESOLIS OF THE BARD TIBE NINOOWS 
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CPU times in seconds on IBM 3033 



2. Soft lime Windows 



All of the methods tested for soft time windows were 
able to find some answer to every problem within reasonable 
computing time, except for two instances with the exact 
algorithm. With the Nearest Neighbor heuristic, the quality 
of solution is not desirable. In general, the lower time 
window percentage proi/lems have lower solution quality. As 
in the hard window case, on test problem [4], the results of 
the Necirest Neighbor heuristic coincide with the uncon- 
strained ISP heuristic, because this problem itself was 
constructed by a Nearest Neighbor. 

As in the hard time window problem, SCCO and SCAO 
generally find an optimal solution except for one problem 
with 50 percent time windows. In test problem [1] with 50 
percent time windows, the SCCO and SCAO values were 215.686, 
165.544 respectively. The exact algorithm could not solve 
the two test problems with 50 percent time windows within 
180 seconds. The reason is that the solutions of the state 
space relaxations have many subtours and it takes a long 
time to eliminate these sub tours. 

With both hard and soft time windows, the results 
are sensitive tc the percentage, width and position of the 
time windows. In most problems, the fewer time windows 
there are, the lower the accuracy of the heuristics. 
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COfiPOTATIOMAL RESOLIS OF THE SOFT TIME WINDOWS 
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VI. COHCLOSIO MS AMD B ECOMM iHD A TlO HS 



This thesis has presented some heuristics and exact 
algorithms for the solution of traveling salesman problem 
with tine window constraints. We considered two different 
kinds of time window constraints ; hard time windows and 
soft time windows. Hard time windows are inviolable, 
whereas soft windows may be violated at a cost. 

For both hard time windows and soft time windows, we 
developed some new heuristics, SCCC and SCAO, which are 
modifications of Stewart’s unconstrained TSP heuristics 
[Ref. 6] CCCO and CCAO. Also for the hard time window only, 
we developed the SLACK heuristic. We also developed an 
exact algorithm for both hard and soft window using state 
space relaxation dynamic programming and branch and bound as 
proposed by Christofides et al. [Ref. .5]. 

The procedures were shown to be effective on some moder- 
ately small sized problems. A Nearest Neighbor heuristic 
was also developed, but it was often unable to solve the 
problem with hard time windows, and it found very low 
guaiity solutions with soft time windows. Tnis experience 
is consistent with tne findings of ethers [Ref. 7] who 
determined that the Nearest Neighbor heuristic does not 
perform well cn the unconstrained TSP. 

The SCCO and SCAO are generally effective cn most of the 
small sized problems we tested, except for the problems in 
which less than half the nodes have time windows. Further 
research is needed in order to satisfactorily solve these 
problems. Another problem difficulty that may reguire more 
research is dealing with wider time windows. 

The SLACK heuristic which is used only with hard time 
windows is slightly slower than the other heuristics. 
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Particularly, in the lower time window percentage problems, 
the accuracy becomes lower. 

The exact algorithm succeeded in solving 14 of the 16 
test problems to optimality, but it was too slow to use in 
most of the lower time window percentage problems. This 
algorithm's performance also depends upon the quality of the 
upper bound which is obtained from the heuristic. 
Additional research is needed to reduce computation time, 
but • a working program for at least some problems has 
resulted from this effort. 
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MEESDIX A 
TEST PBOBLEH £ 1 ] 
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Depot cc-ordinates : (12,10) 

problem source : Sedgewick [Be£. 19] 
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MPEHDIX b 



lESI PfiOBLEH [2] 



node 


X 


y 


1 no de 


X 


y 


1 


2 95 


272 


1 2 


26 7 


242 


2 


301 


258 


13 


259 


265 


3 


309 


260 


14 


315 


233 


4 


217 


274 


1 5 


329 


252 


5 


218 


267 


16 


318 


252 


6 


282 


267 


17 


329 


224 


7 


242 


249 


18 


2b7 


213 


8 


230 


262 


1 9 


275 


192 


9 


249 


268 


20 


303 


201 


10 


256 


267 


2 1 


208 


217 


1 1 


265 


257 


22 


326 


181 



Depot co-ordinates : (326, 1 81) 

problem source : Ellon et al- [Ref. 21] 
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APPEHDIX C 
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£3] 
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source : Ellon et 


al. [fief. 


21 ]. 
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APPEBCI X D 
TEST PfiOBLEM £5] 



I node X y | node x y | node x y | node x y | 



1 


37 


52 


2 


49 


49 


3 


52 


64 


4 


20 


26 


5 


40 


30 


6 


21 


47 


7 


17 


63 


8 


31 


6 2 


9 


52 


33 


10 


51 


21 


11 


42 


4 1 


12 


31 


32 


13 


5 


25 



14 


12 


42 


15 


36 


16 


16 


52 


41 


17 


27 


23 


18 


17 


33 


19 


13 


13 


20 


57 


58 


21 


62 


42 


22 


42 


57 


23 


16 


57 


24 


8 


52 


25 


7 


38 


26 


27 


68 



27 


30 


48 


28 


43 


67 


29 


58 


48 


30 


58 


27 


31 


37 


69 


32 


38 


46 


33 


46 


10 


34 


6 1 


33 


35 


62 


63 


36 


63 


69 


37 


32 


22 


38 


45 


35 


39 


59 


15 



40 


5 


6 


41 


10 


17 


42 


21 


10 


43 


5 


64 


44 


30 


15 


45 


30 


10 


46 


32 


39 


47 


25 


32 


48 


25 


55 


49 


48 


28 


50 


56 


37 



Depot co-ordinates : (30,40) 

problem source : Eilon et al. £Ref .21]. 
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A PPEaJl X E 
TEST PBOBLEH £6] 



1 


node 


X 


y 


1 node 


X 


y 


1 node 


X 


y 


1 


node 


r X 


y 


1 


1 


1 


22 


22 


20 


66 


14 


39 


30 


60 


1 


58 


40 


60 


1 


1 


2 


36 


26 


21 


44 


13 


40 


30 


50 


1 


59 


70 


64 


1 


1 


3 


21 


45 


22 


26 


13 


41 


12 


17 


1 


60 


64 


4 


1 


1 


4 


45 


35 


23 


11 


28 


42 


15 


14 


1 


6 1 


36 


6 


1 


1 


5 


55 


20 


24 


7 


43 


43 


1 6 


19 


1 


62 


30 


20 


i 


1 


6 


33 


34 


25 


17 


64 


44 


2 1 


48 


1 


63 


20 


30 


1 


1 


7 


50 


50 


26 


41 


46 


45 


5C 


30 


i 


64 


15 


5 


1 


I 


8 


55 


45 


27 


55 


34 


46 


5 1 


42 


1 


65 


50 


70 


1 


1 


9 


26 


59 


28 


35 


16 


47 


5C 


15 


1 


66 


57 


72 


1 


1 


10 


40 


66 


29 


52 


26 


48 


48 


21 


1 


67 


45 


42 


1 


1 


11 


55 


65 


30 


43 


26 


49 


1 2 


38 


1 


68 


38 


33 . 


1 


1 


12 


35 


51 


31 


31 


76 


50 


1 5 


56 


1 


69 


50 


4 


1 


1 


13 


62 


35 


32 


22 


53 


51 


29 


39 


1 


70 


66 


3 


1 


1 


14 


62 


57 


33 


26 


29 


52 


54 


38 


1 


71 


59 


5 


1 


1 


15 


62 


34 


34 


50 


40 


53 


55 


57 


1 


72 


35 


6 0 


1 


1 


16 


21 


36 


35 


55 


50 


54 


67 


4 1 


1 


73 


27 


24 


1 


1 


17 


33 


44 


3b 


54 


10 


55 


10 


70 


1 


74 


40 


20 


1 


1 


18 


9 


56 


37 


60 


15 


56 


6 


25 


1 


75 


40 


37 


1 


1 


19 


62 


48 


38 


47 


66 


57 


65 


27 


1 








1 



Depot co-ordiiiates ; (40,40) 

problem source ; Ellon et al.-[Bef -21]. 
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IPPEMDIX F 



TEST PBOBLEH FOE THE SCCO 



1 

1 


node 


X 


y 


time 

l(i) 


window 

u(i) 


1 node 

1 


X 


y 


time 

l(i) 


window 

u(i) 


I 

I 


1 


1 


3 


9 


— 


- 


1 11 


10 


13 


10 


17 


I 


1 


2 


1 1 


1 


- 


- 


1 12 


16 


1 4 


2 


9 


1 


1 


3 


6 


8 


27 


36 


1 13 


15 


2 


- 


- 


1 


1 


4 


4 


3 


37 


45 


1 14 


13 


16 


5 


13 


1 


1 


5 


5 


15 


- 


- 


I 15 


2 


12 


- 


- 


1 


1 


6 


8 


11 


12 


23 


1 16 


12 


1 0 


- 


- 


1 


1 


7 


6 


4 


35 


43 


1 










1 


J 


8 


7 


4 


42 


49 


1 










1 


1 


9 


9 


7 


58 


68 












1 


1 


10 


14 


5 


53 


64 


1 










1 



Depot co-ordinates : (12,10) 

problem source 

node locations : Sedgewick [Ref. 19] 

time windows : see Chapter V. 
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APP ESDI! G 
TEST PH 0 BIBS C ] 



1 

J 


node 


X 


y 


time window 
l(i) u(i) 


1 node 

1 


X 


y 


time 

l(i) 


1 


1 


3 


9 


25 


32 


( 11 


10 


13 


10 


1 


2 


11 


1 


46 


53 


I 12 


16 


14 


2 


1 


3 


6 


8 


27 


36 


1 13 


15 


2 


51 


1 


4 


4 


3 


37 


45 


i 14 


13 


16 


5 


1 


5 


5 


15 


18 


28 


1 15 


2 


12 


22 


J 


6 


8 


11 


14 


23 


i 16 


12 


10 


- 


1 


7 


6 


4 


35 


43 


1 








1 


8 


7 


4 


42 


49 


1 








1 


9 


9 


7 


58 


68 


1 








1 


10 


14 


5 


53 


64 


1 









Depot co-ordinates : (12,10) 

CL = 2. 0, CU =5.0 



window 
u (i) 

17 

9 

59 

13 

50 



problem source 

node locations : Sedgewick [Ref. 19] 

tine windows : see Chapter V. 



85 



APPE8DII H 



TEST PEOBLEM C 1“2 ] 



1 


node 


X 


y 


time 


window 


1 node 


X 


y 


time 


window 


1 


1 








l(i) 


u{i) 


1 






l(i) 


u(i) 


1 


1 


1 


3 


9 


25 


32 


1 11 


10 


13 


1 0 


17 


1 


1 


2 


11 


1 


46 


53 


1 12 


16 


14 


2 


9 


1 


1 


3 


6 


8 


27 


36 


1 13 


15 


2 


51 


59 


1 


1 


4 


4 


3 


- 


- 


1 14 


13 


16 


5 


13 


1 


J 


5 


5 


15 


1 8 


28 


1 1 5 


2 


12 


22 


30 


1 


i 


6 


8 


11 


14 


23 


1 16 


12 


10 


- 


- 


1 


1 


7 


6 


4 


35 


43 


1 










1 


1 


8 


7 


4 


42 


49 


1 










1 


1 


9 


9 


7 


58 


68 


1 










1 


1 


10 


1 4 


5 


53 


64 


J 










1 



Depot co-ordinates ; (12,10) 

CL = 2.0, CO = 5.0 
problem source 

node locations ; Sedgewick [Bef. 19] 
time windows : see Chapter V. 
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APPEIDIX I 
TEST PBOBLEH [1-3] 



1 

1 


node 


X 


y 


time 

l(i) 


window 

u(i) 


1 node 

1 


X 


y 


time 

l(i) 


window 
u (i) 


1 

J 


1 


1 


3 


9 


— 


— 


1 11 


10 


13 


1 0 


17 


1 


1 


2 


1 1 


1 


46 


53 


1 12 


16 


14 


2 


9 


1 


1 


3 


6 


8 


27 


36 


1 13 


15 


2 


51 


59 


I 


1 


4 


4 


3 


37 


45 


1 14 


13 


16 


5 


13 


1 


1 


5 


5 


15 


18 


28 


1 15 


2 


12 


22 


30 


( 


1 


6 


8 


11 


14 


23 


1 16 


12 


10 


- 


- 


1 


1 


7 


6 


4 


35 


43 


1 










1 


1 


8 


7 


4 


42 


49 


i 










1 


1 


9 


9 


7 


- 


- 


1 










1 


1 


10 


14 


5 


- 


- 


1 










I 



Depot co-ordinates : (12,10) 

CL = 2. 0, CU =5.0 
problem source 

node locations : Sedgewick [Ref. 19] 

time windows : see Chapter V. 



87 



APPEBDI X J 
TEST PBOBXEM [ 1-4 ] 



1 

1 


node 


X 


Y 


time window 
l[i) u(i) 


i node 

J 


X 


y 


time window 
l(i) a{i) 


i 

1 


1 


1 


3 


9 


25 


32 


1 1 1 


10 


13 


10 


17 


1 


1 


2 


1 1 


1 


46 


53 


1 12 


16 


14 


- 


- 


1 


1 


3 


6 


8 


- 


- 


1 13 


15 


2 


51 


59 


1 


1 


4 


4 


3 


37 


45 


1 14 


13 


16 


5 


13 


1 


1 


5 


5 


15 


- 


- 


1 15 


2 


12 


- 


- 


1 


1 


6 


8 


11 


14 


23 


1 16 


12 


10 


- 


- 


i 


1 


7 


6 


4 


- 


- 


1 










1 


1 


8 


7 


4 


- 


- 


1 










1 


1 


9 


9 


. 7 


58 


68 


i 










1 


1 


10 


14 


5 


— 


- 


1 










1 



Depot Co-ordicates : {12,10) 

CL = 2. 0, CU =5.0 
Problem Source 

node locations : Sedgewick [Ref. 19] 
time windows : see Chapter V. 
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APPE HDI X K 
TEST PHOBLEM [2-1] 



i 

1 


node 


X 


y time window 
l(i) u(i) 


1 

1 


node 


X 


y 


time window 
l(i) u(i) 


1 


1 


295 


27 2 


125 


135 


1 


12 


267 


242 


170 


179 


1 


2 


301 


25 8 


110 


1 18 


1 


13 


259 


265 


193 


202 


1 


3 


309 


260 


102 


1 10 


J 


14 


315 


233 


57 


67 


1 


4 


217 


274 


242 


250 


1 


15 


3 29 


252 


81 


89 


1 


5 


218 


27 8 


23 9 


246 


1 


16 


3 18 


252 


90 


98 


1 


6 


282 


26 7 


141 


149 


1 


17 


3 29 


224 


40 


49 


1 


7 


242 


24S 


279 


2 86 


1 


18 


267 


213 


382 


393 


1 


8 


230 


26 2 


261 


271 


1 


19 


275 


192 


404 


413 


1 


9 


249 


26 8 


206 


2 15 


1 


20 


303 


201 


432 


442 


1 


10 


256 


26 7 


200 


208 


1 


21 


208 


217 


3 23 


332 


1 


11 


265 


25 7 


183 


193 


1 


22 


3 26 


181 










Depot co-ordinates : 


(326, 181) 












CL = 


2.0, 


CU = 


5.0 















problem source 

node locations : Eilon et al- [Ref. 21], 
time windows ; see Chapter V. 
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APPE HDIX L 
TEST PROBLEM [2-2] 



J 

1 


node 


X 


y time window 
l(i) u(i) 


1 node 


X 


y 


time window 
l(i) u(i) 


1 

i 


1 


1 


295 


272 


125 


135 


1 12 


267 


242 


170 


179 


J 


1 


2 


301 


25 8 


1 10 


1 18 


1 13 


259 


265 


- 


- 


1 


1 


3 


309 


26 0 


102 


1 10 


1 14 


315 


233 


57 


67 


1 


1 


4 


217 


274 


242 


250 


1 15 


3 29 


2 52 


81 


89 


1 


1 


5 


218 


27 8 


239 


246 


1 16 


3 18 


2 52 


90 


98 


1 


1 


6 


282 


267 


141 


1 49 


1 17 


329 


224 


40 


49 


I 


1 


7 


242 


249 


279 


286 


1 18 


2 67 


213 


382 


393 


1 


I 


8 


230 


262 


261 


271 


1 19 


275 


192 


404 


413 


i 


1 


9 


249 


26 8 


206 


2 15 


1 20 


3 03 


201 


432 


442 


i 


1 


10 


256 


26 7 


200 


208 


1 21 


208 


217 


- 


- 


1 


1 


11 


265 


257 


183 


193 


1 22 


326 


1 81 


- 


- 


1 



Depot co-ordinates : (326, 181) 

CL = 2. 0, CU =5.0 
problem source 

node locations ; Eilon et al- [Ref. 21] 
time windows : see Chapter V . 
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TEST PHOBIEH £2-3] 



1 


node 


X 


y 


time 


window 


1 


node 


X 


y 


time 


windo w 


1 


1 










u(i) 


1 








1 (i) 


u (i) 


J 


1 


1 


295 


272 


125 


135 


1 


12 


267 


242 


170 


179 


] 


1 


2 


301 


25 3 


- 


- 


1 


13 


2 59 


265 


- 


- 


1 


1 


3 


309 


26 0 


- 


- 


1 


14 


315 


233 


57 


67 


1 


1 


4 


217 


274 


242 


2 50 


1 


15 


3 29 


252 


81 


89 


1 


1 


5 


218 


27 8 


239 


246 


1 


16 


3 18 


252 


90 


95 


i 


1 


6 


282 


26 7 


141 


149 


1 


17 


3 29 


224 


40 


49 


1 


1 


7 


242 


24 9 


- 


- 


1 


18 


267 


213 


3 82 


393 


1 


1 


8 


230 


262 


- 


- 


1 


19 


275 


192 


- 


- 


1 


1 


9 


249 


26 8 


206 


215 


1 


20 


3 03 


201 


432 


442 


I 


1 


10 


256 


26 7 


20 0 


208 


1 


21 


2 08 


217 


3 23 


332 


J 


1 


11 


265 


257 


183 


193 


1 


22 


326 


181 


- 


— 


1 



Depot co-ordinates : (326, 181) 

CL = 2. 0, CU = 5. 0 
problem source 

node locations : Zilon et al. [Ref. 21]. 
tine windows : see Chapter V. 
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APP EHDIX N 
TEST PB0B1E8 [2-4] 

























node 


X 


y 


time 


window 


1 node 


X 


y 


time window 


i 










u(i) 








1 (i) 


u (i) 


1 


1 


295 


272 


- 


- 


1 12 


267 


2 42 


170 


179 


1 


2 


301 


25 8 


110 


1 18 


1 13 


259 


265 


- 


- 


j 


3 


309 


26 0 


- 


- 


1 14 


315 


233 


57 


67 


j 


4 


217 


274 


242 


2 50 


1 15 


329 


252 


- 


- 


1 


5 


218 


27 8 


239 


246 


1 16 


318 


252 


90 


98 


1 


6 


282 


26 7 


141 


149 


1 17 


329 


224 


- 


- 


1 


7 


242 


24 9 


279 


286 


1 18 


267 


213 


- 


- 


1 


8 


23 0 


26 2 


- 


- 


1 19 


275 


192 


404 


413 


1 


9 


249 


268 


206 


215 


1 20 


3 03 


201 


432 


442 


1 


10 


256 


26 7 


- 


- 


1 21 


208 


217 


- 


- 


1 


11 


265 


257 


- 


- 


1 22 


3 26 


181 


- 


- 


1 



Depot cc-orainates : (326,181) 

CL = 2. 0, CU =5.0 
problem source 

node locations : Lilon et al. [Ref. 21]. 
tine windows : see Chapter V. 
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APP EJDIX 0 
TEST PSOBLEM [3-1] 



J 

1 


node 


X 


y time window 
l(i) u(i) 


1 node 


X 


y 


time window 
1(1) u(i) 


1 

1 


1 


1 


151 


264 


196 


204 


1 12 


156 


217 


105 


118 


i 


1 


2 


159 


261 


185 


193 


1 13 


129 


214 


259 


271 


1 


1 


3 


130 


254 


217 


2 25 


1 14 


1 46 


2 08 


2 


10 


1 


1 


4 


128 


252 


222 


234 


1 15 


164 


208 


92 


105 


1 


1 


5 


163 


247 


174 


185 


1 16 


141 


206 


10 


19 


1 


1 


6 


146 


24 6 


142 


154 


1 17 


1 47 


193 


54 


68 


1 


( 


7 


161 


242 


166 


173 


1 18 


1 64 


193 


79 


89 


I 


1 


8 


142 


239 


13 1 


142 


1 19 


1 29 


139 


30 


38 


1 


1 


9 


163 


236 


159 


165 


1 20 


155 


185 


67 


75 


1 


1 


10 


148 


232 


123 


131 


1 21 


139 


182 


40 


53 


1 


1 


11 


128 


23 1 


242 


2 53 


1 22 


1 45 


215 


— 


- 


1 



Depot co-ordinates ; (145,2 15) 

CL = 2. 0 , CU = 5.0 
problem source 

node locations ; Eilon et al. [Eef. 21], 
time windows : see Chapter V. 
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APPE 8 D1X P 
TEST PBOBLEM [3-2] 



1 

1 


node 


X 


y time window 
l(i) u(i) 


I node 


X 


y 


time window 
1 (i) u(i) 


1 

j 


1 


1 


151 


264 


196 


204 


1 12 


156 


217 


105 


118 


1 


1 


2 


159 


261 


185 


193 


1 13 


1 29 


214 


- 


- 


1 


1 


3 


130 


254 


217 


225 


1 14 


1 46 


2 08 


2 


10 


1 


1 


4 


128 


25 2 


222 


234 


I 15 


164 


208 


92 


105 


1 


1 


5 


163 


247 


174 


185 


I 16 


141 


206 


10 


19 


I 


1 


6 


146 


246 


142 


154 


1 17 


1 47 


193 


54 


68 


1 


1 


7 


161 


242 


166 


173 


1 13 


164 


1 93 


79 


89 


1 


1 


8 


142 


23 9 


131 


142 


1 19 


129 


189 


30 


38 


1 


1 


9 


163 


236 


159 


1 65 


1 20 


155 


185 


67 


75 


1 


1 


10 


148 


23 2 


123 


131 


1 21 


139 


182 


- 


- 


1 


1 


11 


128 


23 1 


242 


253 


1 22 


145 


215 


- 


- 


J 



Depot co-ordinates : (145,215) 

CL = 2- 0 , CU = 5.0 
problem source 

node locations : Eilon et al [ Eef . 21]- 

titte windows : see Chapter V- 
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APPEB DIX Q 
TEST PROBLEM £3-3] 



1 

1 


node 


X 


y 


time window 
l(i) u(i) 


1 node 


X 


y 


time window 
1 (i) u(i) 


1 

1 


1 


1 


151 


264 


196 


204 


1 12 


1 56 


217 


105 


1 18 


1 


1 


2 


159 


261 


- 


- 


1 13 


129 


214 


- 


- 


1 


1 


3 


130 


254 


- 


- 


1 14 


1 46 


2 08 


2 


10 


1 


1 


4 


128 


252 


222 


2 34 


1 15 


164 


208 


92 


105 


1 


1 


5 


163 


24 7 


174 


185 


1 16 


141 


2 06 


10 


19 


i 


1 


6 


146 


24 6 


142 


154 


1 17 


147 


193 


54 


68 


1 


1 


7 


161 


24 2 


- 


- 


1 18 


164 


1 93 


79 


89 


i 


1 


8 


142 


23 9 


- 


- 


i 19 


1 29 


189 


- 


- 


1 


1 


9 


163 


236 


159 


165 


1 20 


1 55 


185 


67 


75 


1 


1 


10 


148 


23 2 


123 


131 


1 21 


139 


182 


40 


53 


1 


1 


11 


128 


23 1 


242 


2 53 


1 22 


1 45 


215 


- 


- 


1 



Depot co-ordinates : (145,2 15) 

CL = 2. 0 , CU = 5.0 
problem source 

node locations ; Eilon et al. [Ref. 21]. 
time windows : see Chapter V. 
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APPE BDIX fi 
TEST PEOBLEH £3-4] 



1 

1 


node 


X 


y 


time window 
l(i) u(i) 


J node 


X 


y 


time window 
l(i) u{i) 


J 

1 


J 


1 


151 


264 


- 


— 


1 12 


156 


217 


105 


113 


1 


1 


2 


159 


26 1 


. 185 


193 


1 13 


129 


214 


- 


- 


1 


1 


3 


130 
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Depot co-ordinates : {145,215) 

CL = 2. 0 , CU = 5.0 
problem source 

node locations : Eilon et al. [Ref. 21]. 
tine windows : see Chapter V. 
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Depot co-ordinates ; (145,2 15) 

CL = 2.0, CD = 5.0 
problem source 

node locations : lilon et al [Eef. 21]. 
time winaows ; see Chapter V. 
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1££EJDIX T 
TEST PROBIEH 
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Depot co-ordinates : (145,215) 

CL = 2. 0, CU = 5.0 
problem source 

node locations : Eilon et al. [3ef. 21]. 
time windows : see Chapter V. 
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APPEMDI X 0 
TEST PBOBLEH [4-3] 
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Depot cc-ordinates : (145,2 15) 

CL = 2.0, CU = 5.0 
problem source 

node locations ; Eilon et al- [Ref. 21] 
time windows : see Chapter V. 
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APPEHDIX ¥ 



TEST PHOBLEH £4-4] 
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Depot co-ordinates : (145,215) 

CL = 2. 0, CU = 5 .0 
problem source 

node locations : Eilon et al. [Ref- 21]. 
time windows : see Chapter V. 
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